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I. INTRODUCTION 


This Final Report covers the period from August 1993 to December 1996. 

During this time three Master's degrees were awarded to the students who 
were supported by the NASA Grant. A fourth student who was supported 
partially by this Grant will be submitting his Master's Thesis by April 1997. 


A visiting Professor from China was awarded a Post Doctoral Fellowship 
under this Grant for two years. One of the Master's degree holders has 
completed all course requirements for his Ph.D. and is working toward his 
Dissertation. This dissertation is expected to be submitted to the Graduate 
School at Howard University by March 1998. This student was also given 
financial help under this Grant. 


Five papers were sent for publication. One of these was published in the 
IEEE - MTT's Proceedings and three were published in the International 
Journal of Electronics. The fifth one was sent to reviewers for evaluation. 


The lists of Master's Theses, their headings and their Abstracts are given in 
the following pages. Also are given the lists of Publications and their copies. 




II. List of Master's Theses 


1. "Coupling single-mode fiber to uniform an dtapered thin-film Waveguide" 

Gadi Jagannath, December 1994. 

2. "A weakly guiding approximation for the propagation characteristic of single 

and double-clad optical fibers". 

Trevor Correia, August 1994. 

3. "Rectangular patch microstrip antenna - analysis and design" 

James Wesley Hall, July 1994. 


A fourth Master's Thesis will be submitted in April 1997. 


A Ph.D. Dissertation will be submitted in March 1998. 




III. Headings and copies of abstracts of the Master's theses. 




HOWARD UNIVERSITY 


Coupling Single-mode Fiber to Uniform 
and Tapered Thin-film waveguide. 


* 


A Thesis Submitted to the Faculty of the 
Graduate School of Arts and Sciences 


of 


HOWARD UNIVERSITY 
partial fulfillment of the requirement for the degree of 


Master of Engineering 
by 

Jagannath V. Gadi 


Washington, D.C. 
December, 1994. 



ABSTRACT 


Theoretical derivations, computer analysis and test data are provided to demonstrate 
that the cavity model is a feasible one to analyze thin-substrate, rectangular-patch microstrip 
antennas. Seven separate antennas were tested 

Most of the antennas were designed to resonate at L-band frequencies (1-2 GHz). 
One antenna was designed to resonate at an S-band (2-4 GHz) frequency of 2.025 GHz 
All dielectric substrates were made of Duroid, and were of varying thicknesses and relative 
dielectric constant values. 

Theoretical derivations to calculate radiated free space electromagnetic fields and 
antenna input impedance were performed. MATHEMATICA 2.2 software was used to 
generate Smith Chart input impedance plots, normalized relative power radiation plots and 
to perform other numerical manipulations. Network Analyzer tests were used to verify the 
data from the computer programming (such as input impedance and VSWR) Finally, tests 
were performed in an anechoic chamber to measure receive-mode polar power patterns in the 
E and H planes. 

Agreement between computer analysis and test data is presented The antenna with 
the thickest substrate (e r = 2.33 , 62 mils thick) showed the worst match to theoretical 
impedance data This is anticipated due to the fact that the cavity model generally loses 
accuracy when the dielectric substrate thickness exceeds 5% of the antenna s free space 
wavelength [ 1,2,3] A method of reducing computer execution time for impedance 
calculations is also presented. 
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HOWARD UNIVERSITY 


WEAKLY GUIDING APPROXIMATION FOR THE 
PROPAGATION CHARACTERISTICS OF 
SINGLE AND DOUBLE-CLAD OPTICAL FIBERS 


A Thesis 

Submitted to the Faculty of the 
Graduate School of Arts and Sciences 

of 

HOWARD UNIVERSITY 


in partial fulfillment of 
the requirements for the 
degree of 


MASTER OF ENGINEERING 
Department of Electrical Engineering 


by 

Trevor C. Correia 


* 


Washington, D. C. 
August 1994 



ABSTRACT 


Recent advances in optical fiber technology have made it possible to develop optical 
communication systems with very large bandwidths and very low loss. Optical systems are 
routinely tested and installed with repeaterless spacings of several hundred kilometers, a fact 
that attests to the very low loss inherent to the transmission media. Although there is a 
manufacturing capability for producing high performance optical fibers, there is no concise 
method of accurately predicting their propagation characteristics. For the simplest form of 
fiber, the standard single-clad step-index profile, relatively simple analytical procedures give 
accurate results. However, for the high performance modified structures, exact analytical 
solutions do not exist so their analysis relies on approximate solutions. Because of the rigors 
associated with accurately analyzing optical fibers with modified characteristics, new 
mathematical procedures are constantly being developed to simplify the analysis, and, when 
all else fails, the procedures turn to empirical methods. 

In this Thesis, an approximation method, known as the weakly guiding approximation, 
is utilized to find the propagation characteristics of single and double-clad optical fibers. 
Using a combination of numerical and graphical solutions, data are provided to compare the 
results from this method with those achieved by other methods. It is shown here that, without 
the complexities associated with other procedures, this method gives fairly accurate results. 



HOWARD UNIVERSITY 


RECTANGULAR-PATCH MICROSTRIP 
ANTENNA ANALYSIS AND DESIGN 


* 


A Thesis 

Submitted to the Faculty of the 
Graduate School of Arts and Sciences 


of 


HOWARD UNIVERSITY 


in partial fulfillment of 
the requirements for the 
degree of 


MASTER OF ENGINEERING 
Department of Electrical Engineering 
by 

James Wesley Hall 


* 


Washington, D.C. 
August 1994 



ABSTRACT 


Theoretical derivations, computer analysis and test data are provided to demonstrate 
that the cavity model is a feasible one to analyze thin-substrate, rectangular-patch microstrip 
antennas. Seven separate antennas were tested . 

Most of the antennas were designed to resonate at L-band frequencies (1-2 GHz). 
One antenna was designed to resonate at an S-band (2-4 GHz) frequency of 2.025 GHz. 
All dielectric substrates were made of Duroid, and were of varying thicknesses and relative 
dielectric constant values. 

Theoretical derivations to calculate radiated free space electromagnetic fields and 
antenna input impedance were performed. MATHEMATICA 2.2 software was used to 
generate Smith Chart input impedance plots, normalized relative power radiation plots and 
to perform other numerical manipulations. Network Analyzer tests were used to verify the 
data from the computer programming (such as input impedance and VSWR) . Finally, tests 
were performed in an anechoic chamber to measure receive-mode polar power patterns in the 
E and H planes. 

Agreement between computer analysis and test data is presented . The antenna with 
the thickest substrate (e r = 2.33 , 62 mils thick) showed the worst match to theoretical 
impedance data. This is anticipated due to the fact that the cavity model generally loses 
accuracy when the dielectric substrate thickness exceeds 5% of the antenna's free space 
wavelength [ 1,2,3], A method of reducing computer execution time for impedance 
calculations is also presented. 
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IV. List of Publications 


1. "Coupling Single-Mode Fiber to Uniform and Symmetrically 

Tapered Thin-Film Waveguide Structures using Gadolinium 
Gallium Garnet". Gadi, Yalamanchili and Shahid, SBMO/ 

IEEE MTT-s 1995 Proceedings.pp 926-931. 

2 . "Rectangular Waveguides With Two Conventional and Two 

Superconducting Walls". Yalamanchili,Qiu and Wang. 
International Journal of Electronics, 1995. 

3. "Scattering by the Transition Junction Between Conventional 

and a High-Temperature Superconducting Waveguide". 
Yalamanchili, Qiu and Wang. International Journal of Electronics, 
Vol. 76, No. 6, pp 1 173 - 1180, 1994. 

4 . "Meissner Model of Superconducting Rectangular Waveguides". 

Wang, Qiu and Yalamanchili. International Journal of Electronics, 
Vol. 76, No. 6, pp 1151 - 1171, 1994. 

*5. “Single Index Summation Procedure for Calculating Antenna Input 
Impedance”. Yalamanchili and Hall. 1996. 


* This paper has been sent to reviewers for evaluation. 




V. Copies of Published Papers 




COUPLING SINGLE-MODE FIBER TO UNIFORM AND SYMMETRICALLY 

tapemdthVn-film waveguide structures using gadolinium 

OALLIUM garnet 


Jagannath Gadi, Raj Yalamanchili and Mohammad Shahid 
Howard University 
Department of Electrical Engineering 
2300, 6 th Street, N W 
Washington, DC 20059 


Abstract 

The need for h,gh efficiency components has grown significantly due to the expanding role of fiber opi,c 
communications for various applications Imegraied optics is in a slate of metamorphosis and there are many 
problems awaiUng solutions One of the main problems being the lack of a simple and efficient method of 
smgie-mode fibers io thin-film devices for integrated optics In this paper, optical coupling between a 
single-mode fiber and a uniform and tapered thin-film waveguide is theoretically modeled and analyzed A novel 
tapered structure presented in this paper ,s shown to produce perfect match for power transfer 

1. INTRODUCTION 

Integrated Optics has come a long way since early I970's. II is poised al the threshold of making a big 
impact in everyday world Integrated optics has drawn into several disciplines, such as computers and 
microwave integrated circuit technology resulting in new fabncauon technologies The field is now in a srate of 
flux and there are still many problems awaiting solutions One of ffie main problems being die lack of a simple 
and efficient method of coupling opUcal fiber 10 Ihin-film devices fo, imegrated opucs Although a number of 
connectors are available in the market, the efficiency and reproducibility are low In this paper, a novel uniform 
and symmetrically tapered structure are analyzed mathematically and ihe resells presented 

There are certain papers published by resrarchers in the area of thin-film couplers that are worthy of note, in that 
they provide the background for the development of the proposed wortt and point towards issues raised b, 
shortcomings of me previous wort: Louisell 1 1 1 invesfigated broadband bi-dirachomrt couplers ,n which he shownd 
mat me phase constants and coupling coefficients vara with disumce along ram coupled transmission lines 
Ulrich 12 | has shown by analysis how ligh. ran be coupled inio a ffii.-f.lm by mrans of a pnsm-film coupler Akira 
Ihaya PI presented a mathemaucal model of a thin-film opucal directional coupler consisting of a three-layered 
deposited glass films on the substrate, with coupling occurring beraiecn me first and the third film Wilson and 
TehHl have shown a mathemaucal modeling of a lapered velocity direcuonal coupler Nelson |5| has theoreurally 
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examined the coupling of single-mode optical waveguides through the use of expanding and contracting tapers 
Juichi Noda et al.(6] have shown a connection between single mode fiber coupler to a Ti diffused LiNo 3 strip 
waveguide The authors Y Cai et al.[7) , have analyzed the coupling characteristics of a uniform structure 

2. THEORETICAL BACKGROUND 

In this paper, we present a novel method of coupling light from a single mode fiber to two structures Uniform 
and symmetrically tapered thin-film couplers. The model presented is an improved version of the uniform five 
layered structure presented by Y Cai et al [7], If a fiber is directly connected to the thin-film. as shown by Y.Cai et 
al 17].. it leads to a large mismatch of the field profiles at the interface and as a result leads to loss of optical 
power. To overcome this optical reflection and radiation losses, they proposed a five layered structure which is 
shown in Figure 1. This consists of a coupling waveguide, buffer layers, thin-film layer and the substrate To avoid 
mismatch of the fiber core and the coupling waveguide field profiles at the interface, the coupling waveguide 
dimensions are designed as 2aX2a. where a' is the radius of the fiber core Gutunann et al(8].. have shown that 
the field distributions of the fundamental mode of the cylindrical fiber core and the garnet material waveguide 
differ slightly if the refractive index differences are smaller Solgcl solution is used at the interface such that any 
slight refractive index mismatch between the fiber core and the coupling waveguide is minimized The buffer layer 
serves the purpose of coupling optical power to the thin-film waveguide of lower refractive index. Gadolinium 
Gallium Garnet material was used as the thin-film material because of the properties exhibited by the material 
The buffer layers and the coupling waveguide are silica doped materials. The required refractive index of buffer 
and coupling waveguides is obtained by doping TiO; with SiO; 

The eigen mode equations of the five layered structure developed ( Figure 1 ) can be obtained by solving with 
w=g=2a. where : w= the height of the coupling waveguide. a= radius of the fiber core. n s =rcfractive index of the 
coupling waveguide ( Garnet material ). n.=refractive index of the buffer layer. n w = refractive index of the thin- 
film material. X=wavelength of the wave ( 1.3um ) and Wl= height of the thin-film coupling waveguide The 
eigenmode equations developed by Y Cai et al[9]..are 

Q a = nlK gx x*n(K gx W/2) ~ (1) 

Q a = V (2* / W(" 2 „ ~ n l)~ K l (2) 

o ' ’ 

* t*n(gK„) 

Q, = JpnTWtf. “ (4) 
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The parameters k^. Q, and k^. Q, are the x and y components of the wavenumbers and represent the 
sinusoidal vanaUon in the x and y directions The propagation constant of the coupling waveguide can be 
determined by solving equations ( 1 to 4 ), 

P\ s = 0} : Mo e 0 ”] s ~K~ K l ( 5 } 

The propagation constant of the thin-film waveguide is obtained by replacing n* with n*. w with W1 and 
with p w in equations! lto5 ) The value of W1 for the thin-film is so chosen such that highest power 
coupling occurs between the coupling waveguide and the thin-film at this value The optimum value W1 is chosen 
such that P w = 3w for no mismatch The thickness of the buffer layer has a direct bearing on the power coupled 
from the coupling waveguide to the thin-film waveguide and an optimum value is chosen such that there is a 
maximum power transfer from the coupling waveguide to the thin-film waveguide. The propagation constant of 
the thin-film vanes because of the change in value of w along the z direction. Therefore the mismatch. M depends 


on the propagation constant difference given by M{z) - 


x _ r gg 


nr 


( 6 ) 


and the maximum fraction of power transferred between the coupling waveguide and the thm-film waveguide with 

a fixed mismatch M. is given by F : =(l+M : )' (7) 

Application of the theory developed by Snyder et al| 10] and Snyder ( 1 1]. to the model analyzed in Figure 1 . gives 

the power coupled into the thin-film for a uniform coupler as 


F =PF : Sin : \{CIF)dz (8) 

0 

where P is the total power introduced. C the coupling coefficient between the center of the coupling waveguide 
and the center of the thin-film waveguide In this paper, we present a novel uniform and symmetrically tapered 
structure ( Figure 2 ) where the modes of a uniform waveguide form a complete set and can propagate 
independently from one another, while the tapered modes are coupled together and adjust their characteristics to 
suit the varying transverse properties of the guiding structure as they are moving along the taper. During the 
process of coupling the light from the uniform coupling waveguide to the thin-film waveguide, the wave is kept in 
the lowest order mode. For the shape shown in this paper.! Winn and Hams [12]).the power coupled at the 

narrow end is given by P 0 = P (l*F*Sin J (C/F^d*) ^ 

0 

In this paper while analyzing the theoretical model of the uniform and tapered coupler we made the assumption 
that the materials used are lossless and therefore there are no Fresnel’s reflections at the interfaces and that there 
are no lossess at the interface of the fiber core and the coupling waveguide. 
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3. RESULTS AND DISCUSSION 


The material used for the analysis of uniform and tapered couplers is Gadolinium Gallium Garnet (GGG) The 
coupling waveguide material used is GGG. Buffer layer material is spun silica which is doped with Ti02 Thin- 
film material is polymerized solgel solution of SiO : and TiO : The refractive index of coupling waveguide n» = 
1.9389 (a ). = 1 3 urn The refractive index of the buffer layers is chosen, as n, = 1.9340 and the refractive index 
of the thin-film material, n, = 1 .9450 Figure 3 shows the plot of C versus d. where C is the coupling coefficient 
and d is the distance from the center of Uie coupling waveguide to the center of the thin-film. This equation given 
by Snyder [11]. is 



The refractive index of GGG versus wavelength is shown in Figure 4 The expression used to calculate the 


refractive index is 



( 11 ) 


where A, and L, are the sellmeier coefficients. given by Wood and Nassau [13] 

By solving the equations ( lto5) with w=2a=g. the propagation constant 3* = 9.3253804 urn 1 was obtained 
Figure 5 shows the propagation constant for the thin-film as a function of the film thickness W 1 The propagation 
constant of the thin-film equals the propagation constant of the coupling waveguide when the thickness 
Wl=0.886|j.m The maximum coupling power is attained at (3* = pw = 9.3253804 pm 1 Figure 6 shows the 
plot of coupled power versus the propagation distance for the uniform coupler of different d s for perfectly 
matching constants (3. = (3w at thickness of thin-film w=0 886 pm. Figure 7 shows the power coupled for 
uniform and tapered couplers with respect to the distance z and the taper starts at length Ll= 700pm. with 
slopes (K) of the symmetrical taper at 0.0001. 0.00009 and 0.00002 Figure 8 shows the effect of the thickness of 
the thin-film and the role played by the taper in the power output. 


4. CONCLUSION 

The novel structure (Figure 2) exhibits broad-bandwidth coupling characteristics and it is easy to manufacture 
since it does not put serious constraints on the accuracy of the coupling length of the coupler. The taper has the 
distinctive advantage of confining the power within the taper such that the output stays approximately near the 
value of the power introduced at the start of the taper The outstanding feature of the taper is the higher efficiency 
as compared to the more conventional devices. 
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Rectangular waveguides with two conventional and two 
superconducting walls 

RAJ YALAMANCHILIf, ZHENG AN QUIf and 
YEN-CHU WANGf 

The propagation properties of TE*"" modes and their dispersion relations in 
rectangular waveguides with two conventional and two superconducting walls, 
derived by using the Meissner boundary conditions on the superconducting walls, 
are presented. In addition to recovering some previously known results, some 
novel results have been obtained: the cut-off wavelength of the dominant TE 10 
mode is greater than that of the conventional TE t0 mode, and the tangential 
electric field and normal magnetic field for the dominant mode TE 10 exist on the 
superconducting surfaces. Expressions for electromagnetic components, surface 
currents, attenuation coefficient, maximum transmitted power, dispersion and 
wave impedance are also presented. 


1. Introduction 

Over several years there has been considerable progress in the development of 
superconducting devices in the microwave and millimetre wave bands, for example, 
low- and high-temperature superconducting microwave filters, resonators, antennas, 
phase shifters, etc. Many superconducting electronic devices have been sought at 
liquid nitrogen temperatures (NisenofT 1988, Van Duzer and Tuner 1981). The 
discovery of high-7^ superconductors (HTS) has fundamentally changed the pro- 
spects of superconductive electronics. The low surface resistance of superconducting 
materials makes possible microwave devices and circuits with very high Q t low insert 
loss and dispersion. Superconducting waveguides have been studied as low T c 
superconductors (< 18 K) (Alaux and Wybouw 1976, Rohner 1978, Zepp et ai 1977, 
Fabre et ai 1981). Experimental and theoretical results have shown that HTS 
waveguides and transmission lines exhibit significantly better performance than that 
of their identical metallic counterparts (Wang et ai 1994, Yalamanchili et ai 1992). 

The propagation properties of transverse electric (TE) modes in high-T c 
(T C >30K) superconductor rectangular waveguides have recently been studied using 
HTS electromagnetic theory and Meissner boundary conditions which are consistent 
with the two-fluid model (Wang et ai 1994). The new boundary conditions give rise 
to different field solutions, hence new waveguide characteristics. This paper describes 
the propagation properties of TE modes in rectangular waveguides with two 
conventional and two superconducting walls (WGCSW) based on the theory 
established by Wang et al. (1994). This kind of waveguide has different character- 
istics from either HTS waveguides or conventional waveguides. The attenuation and 
dispersion are smaller than that of conventional waveguides. The dominant mode in 
the WGCSW is the TE 10 mode, which is the same as conventional waveguides. An 
interesting property is that the cut-off wavelength of the dominant mode is greater 
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rnode'hTlhe WGCSW « « easier ,o exciie ,hc dominan, 
" GCSW ll " W °CSW has two conventional broad walls and 

osupe condueng walls (two simply-connected superconductors) it is easier to 
fabricate than a four-wall superconducting waveguide 

the I "-dirccdL 2 Vrom m .h S , t ' mp ° rtant dispersion Nation of WGCSW is derived along 
wo walls °-o T i W p VC equat,on f and the Meissner boundary conditions on 
surfaacurnsnts^*i nrlnH ****?* for ='«tromagnetic field components and 
normal cur«n V tt? , s “P ercondu «i“« current J u and the surface 

wmw J T b ! ° n two ' flu,d theor y a « given. Various parameters of the 
WGCSW are analysed and compared with conventional waveguides. 

2. Theory 

2.1. Wav e equation and numerical solutions of cut-off wavelength 

This section analyses a WGCSW with two narrow superconducting walls (x=0 
and _v-a) and two broad perfect or good conductor walls {y = Q and vi) as shown 
hoimHa ' A ngorous formulation based on Maxwell’s equations and Meissner 
WGCSW. C ° ° nS ' S t<5 ° bta,n thC " eW electroma 8 net ic field properties in the 

The governing wave equation (Helmholtz equation) and the Meissner boundarv 
WGCSW S a«: '° ng,tUdinal magnelic field H -- of the TE mode in an^ir-filled 


dH. I 

~d^~T L H -- =0 alr=0 ' x = a 


(1) 
(2 a) 


dH 

dn 


“ = 0 aty=0, y = b 


(2 b) 


r> h nUnl h i e p 0PCrat0r ^ '1 thc La P lacian operator in the transverse plane (i e 
u -plane). Propagation in the .-direction is assumed. The parameters a and b are the 


Figure f. 
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wdth and height of the WGCSW, respectively, k c is the cut-off wavenumber, and 
is London penetration depth is for HTS this value varies from |<r*m to KT’m 
( urns 1992). Equation (2a) is the Meissner boundary condition; (2b) is the 
t b he U field ry H COnd ' t ‘ 0n f0 u 'l* convcntionaI wal1 - The important Meissner effect is that 

th!v fi H^ ° n ° Van,S , UP y fr ° m the SUrface of a bulk superconductor, rather 
they decay exponentially into the bulk. The penetration depth X L is the characteristic 

decay length of the magnetic field into a superconductor. It is obvious that when 

V*' thc r°u d term 0f thc Mcissncr houndary condition approaches zero (2a) 
reduces to the boundary condition (*//,/*«) = 0 for the conventional waveguide 

. * f °. re ’ thC , firSt ' Crm of ( 2fl ) 1S related to the properties of the WGCSW Thus 

modd C1SSnCr b ° Undary condltions are relevant to the two-fluid superconducting 

It is assumed that the general solution of the governing equation is 

H z = (C t sin + C 2 sin /rpc) cos k^y ^3^ 

J r ' Cl a " d C * are unknown coefficients; both cannot be equal to zero at the same 

Sidi kf thC o n ° Wn waven umbers along the x and y axes, respectively, 
stitutmg (3) into (2 a) and (2 b) obtains two sets of equations for k x and k,: 


(k* sin k x° c i ~ cos k x a + j c 2 =0 


(4 b) 


2 ~~k x C 2 = 0 


and 


(4 b) 


mn 


b ♦ w "0, 1, 2, 3, . . . 

S m a“ s - ,he *— h “ or lhe -* »U4.) 

-k t cos k^a- 


I 

*2 


i ■ ci n t ,, cos i . sin k^a 

k x sin K x a ; — k. cos k^a _ 


= 0 


( 6 ) 


Therefore the one important dispersion relation in the transcendental equation is 
obtained from thc above determinant as 


tan k z a = 


(Ml) 1 -! 


(7) 


The non-integer roots of this equation describe the dispersion relations along the 

WGCSW A fUnCt h 0 " .° f L r, d ° n de P th and ‘he dimensions of the 

W k G ? 1 V ' A 8r i P r' Ca mCth ° d ,S USed for thc solu ‘*on. as shown in Fig 2 from 
which the cut-off features are obtained. 8 ’ 
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Substituting (4 a) into (3), 

W. = C, ^cos k ir x + -^ sin cos kj 

cos kfX 4- sin k r x) cos k^y 


( 8 ) 


(9) 

It is obvious from (9) that *,*0, otherwise the solution becomes trivial It is worth 
notmg that owing to the non-integral nature, the mode indexes are determined by 
solving (7), which ,s denved from .he Meissner bottndary condittons two 

To avoid risk of confusion with the commonly used TE the sunerscrims nm 

r,U n WGCSW e '[Tl™' Tr" dUC,i " 8 Wa " S * “* 1,!td modes 

" Qn f GCSW . In th , ls case ’ the dominant mode is designated TE 10 (/>= l does not 
mean integer index along the x-axis, instead it is the first root of (7)). The above 
formula, however, does cover conventional waveguides in which H.xcosk tcosJk v 
when /. L -^oo, as is evident from (8). There are no TE 0m modes in 'the WGCSW i£ 
which property it is quite different from conventional waveguides 
Substituting (3) into (1), 


k l=k\ + k 2 y 

Therefore, the cut-off wavelength / c is obtained as 


/, = - 


2n 


( 10 ) 


( 11 ) 


L h h„ nU Tv" Cal T'a ° f Cm ' olr wa,dOT * lhs TE- modes are listed in the 

in “ e WGCSW ?“ ,h ' *avelen e th or the dominant mode 

in the WGCSW (* e % 4 06 mm for a=2mm) is greater than that of TE „ in a 

conventional waveguide with the same dimensions by a factor of 1013 and the 
bandwidth is greater than an octave. 
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2.2. Field components and surface current distributions 

Field component expressions for the TE mode can be determined from the axial 
magnetic field H : by means of the relations between the transverse and axial 
components as follows. 

From (9), let B 0 — kJ. L and H 0 = C 2 for simplicity. Then the axial magnetic field 
component H z becomes 


H Z = H 0 {B 0 cos fcpr + sin k r x) cos k y y (12) 

where H 0 is an unknown coefficient depending on excitations, B 0 is a proportion- 
ality parameter which varies for different modes and penetration depths. The wave 
factor in the form of exp(j cot-yz) is assumed. The real part a of the propagation 
constant y is the attenuation constant, the imaginary part /? is the phase constant and 
oj is the angular frequency. Other field components can be expressed as 


H = 


JPK 


U ( D 


(13a) 




j pk, 

-p- H 0 (B 0 cos kxX + sin k^x) sin k y y 


(13*) 


x ~ k ; ? 


H 0 (B 0 cos k^x + sin k r x) sin 


k,y 


(13c) 




H 0 (B 0 sin k r x — cos k^x) cos k y y 


(13rf) 


where n 0 is the permeability of the material, assumed to be that of free space. 
The wave impedance Z h of the WGCSW is defined as 





Ey_ta/io 
H x /? 


(14) 


Cut-off wavelength (mm) 

Modes 

a- 1 



j = 8 

TE 10 

2 05 (2-0)** 

4-06 (4-0) 

8 07 (8 0) 

16 09 (16 0) 

TE 20 

102 (1-0) 

2-04 (2-0) 

4 07 (4 0) 

8-04 (8-0) 

TE 11 

0*90 (0-89) 

1-79 (1-78) 

3-58 (3-58) 

716 (7-15) 

TE 21 

0-70 (0-70) 

1-42 (1-41) 

2-83 (2-82) 

5 66 (5-65) 

TE J0 

0-67 (0-48) 

1-34 (0-97) 

2-67 ( 1 -94) 

5-34 (3 88) 

TE J1 

0-56 (0-44) 

1-11 (0 89) 

2-22 ( 1 -79) 

4-44 (3-58) 


•For /. L = 10 8 m there is little difference in values. 

••Values in parentheses are the corresponding cut-off wavelengths for the first six TE„ modes in 
conventional waveguides. 


Cut-off wavelengths for some TE pm modes (/ L = 10~ 7 m, a = 2b)*. 
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The instantaneous field expressions for the dominant TE 10 mode in the WGCSW 
can be obtained as: 

H z — H 0 (B 0 cos k r x + sin k r x) cos (wt — fiz) (15a) 

Pk x 

H x = —y H 0 {B 0 sin k^x — cos k^x) sin (c ot — fiz) (15 6) 

* c 


cjfu 0 k x 

E y = , 2 H 0 ( — B 0 sin /c x .x-hcos fc^x) sin (a>f — /fc) (15 c) 

*c 

where, = is the phase constant of the TE 10 mode, and x f is the waveguide 
wavelength of the TE 10 mode. There are three field components in the WGCSW, as 
in conventional waveguides. Note that the superscript ‘10’ for the dominant mode is 
omitted from here on unless otherwise stated. 

From the above discussion, (15 c) demonstrates clearly that, due to B 0 being very 
small, the space distribution of the electric field E y of the TE 10 mode is proportional, 
significantly, to cos k^x along the .x-direction. This means that tangential electric 
fields exist on the superconducting surfaces (* = 0 and .x = a) of the WGCSW. This 
novel property is completely different from that of a perfectly conducting conven- 
tional waveguide since, for H z ^0 and for a finite surface impedance at microwave 
and millimetre wavelengths, the tangential components of E on the superconducting 
walls are finite. The space distribution along the v-direction is the same as for 
conventional waveguides. For details of the field distributions see Fig. 3. 

Equation (156) demonstrates that the space distribution of the magnetic field H x 
is largely proportional to cos k^x, which means that normal magnetic fields exist on 
the surfaces (.x = 0 and x = a) (see Fig. 3). The above results demonstrate that the 
magnetic field lines are not continuous in the WGCSW (Orlando and Delin 1991), 
which is in agreement with the Meissner effect, but the magnetic flux density B is 
continuous in and out of the WGCSW walls because £ = /*<,(//+ A/) (Af is the 
magnetization density in the walls and is not equal to zero). This also is a new 
property in the WGCSW, which does not exist in conventional waveguides. 




Figure 3. (a) Distribution of the electric field £,; (6) distribution of the magnetic field 

for the TE 10 mode. 
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London’s equations and the two-fluid model are used to investigate the current 
distribution on the surfaces (x = 0 and x-a). The total conduction current density 
consists of the superconducting current density J % and the normal electron density / n 
based on the two-fluid model. This model assumes that the conducting electrons in 
the WGCSW walls are divided into two categories. We use a different method to 
describe their distribution. The surface superconducting current /„ due to super- 
electrons is calculated based on the London relation and appropriate boundary 
condition first. Then the normal surface current / so due to normal electrons and the 
losses is calculated using the approximate resistive boundary condition and the 
perturbation approach. 

TL^ superconducting current J % in two isolated walls in which each wall is 
considered as a simply-connected superconductor using a thin-film substrate is given 
as 


J,= ~A ( 16 ) 

A 

This is London’s relation, where A — p 0 Xl, and A (in T m) is the vector potential in 
the WGCSW. Let the vector potential inside the air-filled hollow region in the 
waveguide be A x (subscript a indicates the air-filled WGCSW). This satisfies the 
following equations: 


V + A t = B 

(17) 

k 2 A, + VVA, 

(18) 



jo>/i 0 e 0 


where e 0 is the permittivity of the material. 

Thus, the London relation guarantees that the vector potential is a Teal field’ 
that is completely specified. Equation (18) is called Lorentz’s gauge, which gives rise 
to the relation between the electric field E and the vector potential A % in the air 
region. 

Before solving the current distribution problem, it is important to determine the 
boundary condition at a boundary between the two media (here, the superconduct- 
ing medium and air for walls x = 0 and x=a\ the perfect conductor and air for walls 
y = 0 and y = b). The boundary condition on the tangential components of fields (/, 
and A J for the simply-conneccted superconductor is 

4 ,= -^.. 09 ) 

where the subscript ‘at’ denotes the component tangential to the boundary in the air 
or the applied field region; ‘st’ denotes the component of supercurrent density 
tangential to the boundary in the HTS. Note that J M is the bulk current. Once A t is 
known, the current density /„ and the surface supercurrent density are readily 
calculated. The solution satisfying (17) and (18) is obtained from B as follows: 



8 


R Yalamanchili et al. 


= sin k^x + cos k^x) sin k^y 


(20 a) 


.Mo* 


-4.* = — H 0 (B 0 cos fc^x-sin /c r x) sin 

fi 0 k 

— — f ~ U 0 ( B 0 cos k r x + sin ^.t) sin 


(20 b) 


J/4 

^, y = - ^ y (B 0 sin /r^x + cos fc^x) cos ^ 


/i 0 ** 

+ -p- (B 0 sin — cos /c^x) cos /c y v 


(20 c) 


The surface supercurrent distributions on the isolated walls (.x = 0 and x = a) are 
obtained from (19) as follows: 

(x = 0 plane) 


*^S»t J \%y&j Ai:®; 


A k 2 


H 0 cos k y y — 


j^o k f 
P 


cos 



A 0 a . 

a y J" sin k y>' a z 


(21 a) 


(.x = a plane) 

Asl *^5*>'^y ^ js:^; 

JlW, , ]A 0 k y 1 A 0 X 

A mF" H 0 C ° S ^ J 1 COS ^ fly + /[ Sin (2 1 b) 

where and J uz are the surface supercurrent magnitudes on the walls along the v- 
and ^-directions, respectively. Consider, here, that the superconducting current 
decays exponentially in the walls, so the current is mainly carried by the supercon- 
ductor surface adjacent to the air in the WGCSW. It is easy to see thr r the surface 
supercurrent and bulk current are very large because the amplitude of J s%t is 
proportional to 1// L (A =/i 0 /£). which is of the order of 10 7 . 

All high-7^ materials are type II superconductors. In the Meissner phase, it is 
assumed that the critical field f/ cl ^7*7 to 385 Am -1 (depending on temperature) 
for a HTS sample of YBa 2 Cu 3 O y (88 2 K, axis) (Wu and Sridhar 1990), so if 
//o^I 3 to 660 Am then the theoretical values of supercurrent density are 
estimated to be about lO 6 ^ 8 A cm 2 which is in excellent agreement with other data 
(Burns 1992, Heinen et al. 1991, Miranda et at. 1991, Levenson et al. 1991). The 
distributions of surface supercurrent and normal current on the .x = 0 and .x = a 
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planes of the WGCSW are shown in Fig. 4. (Note that the distributions of the real 
and imaginary parts of the surface supercurrent of the WGCSW are the same.) 

e surface norma! current distribution at its boundaries can be obtained using 
the approximate resistive boundary condition nxH=J tn (Senior 1975) as given by 

(on the ^ = 0 plane) 


An ~ + Anz°z 


- H 0 (B 0 cos kjc + sin k^)a x H 0 (B 0 sin * r r-cos (22 a) 


(on the y = b plane) 

An Anx^x F Anz 

= —H 0 cos kyb(B 0 cos k r x + sin k x x)a x 
j 0k x 

+ ^2 Ho cos k y b(B 0 sin k^x — cos k^xja. 
(on the x = 0 plane) 

m y 


(22 b) 


An - Any a y + Anz<*z = H 0 B 0 cos kja, - J -p^ H 0 B 0 sin k y ya : (22 c) 


(on the x = a plane) 


An Any a y + Anz a z 


— H 0 (B 0 cos k^xi + sin k^) cos k y ya y 

jpk, 

£2 H 0 (B 0 cos k^a — sin k x a) sin k^a. 


(22 d) 


: h r/*"; an ^* n - are the Surfacc n0nnal current ma « nit ud« on the walls alone the 
and z-directions, respect.vely. It is easy to see that the distributions of the surface 
normal currents on the walls .t = 0 and are the same as those of the surfa« 
supercurrents on the walls, which is the result of the two-fluid model. But the 
distributions of current on the y = 0 and y = a surfaces differ from conventional 
WGCSwl 'f 10 different electromagnetic field distributions. The losses of the 
WGCSW come from the surface normal current flowing through the small surface 




Figure 4. Distribution of surface supercurrcnt and surface normal current for the TE‘° 
mode on the walls (a) .r = 0, (6) x=a. 
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impedance. It is predictable that the losses are small because B 0 is very small. The 
distributions of surface current on the WGCSW planes >> = 0 and y = b are shown in 
Fig. 5. 

2.3. Attenuation coefficient 

The attenuation caused by two superconducting walls and two conventional 
walls can be calculated by the perturbation technique as given approximately by 



where f c is the cut-ofT frequency of the TE 10 mode, R, and R n are the surface 
resistances of HTS and metal walls, respectively, and a l is in millimetres. Note that 
the attenuation is not only proportional to because R,oc/ 2 (Romanofsky and 
Bhasin 1992), but also proportional to (R n oc/ 1/2 ). The curve of attenuation 
coefficient ot versus frequency / is shown in Fig. 6 in which R # ^5* Ox 10 28 / 2 Q, 
R n ^2-6x 10~ 7 / l/2 Q (Cu), and a = 26 = 2mm. As far as the attenuation is con- 
cerned, it is greater than that of high-T c superconducting waveguides, but smaller 
than that of conventional waveguides. 

Losses in the WGCSW are generated from three sources: superconductor loss, 
metal loss and dielectric loss. A typical superconducting waveguide used in micro- 
wave and millimetre-wave circuits uses low-loss dielectrics or free space, so the 
dominant loss mechanisms are the superconductor loss and metal loss. 

The superconductor loss is caused by the current of normal electrons in the 
superconducting walls based on the two-fluid model theory. The paired supercon- 
ducting electrons are dissipationless because they cannot be scattered without 
breaking pairs. The microwave surface resistance of superconducting films has been 
found to be as low as a few microohms. These surface resistance values are lower 
than those for copper at the same temperature by a factor of several tens. Such a 
surface resistance gives a small power loss. If the penetration depth x L = 3-6 x 10 7 m 
for the case of conduction along the c-axis and the surface resistance is equal to 
10-3 mfl (60 K) at 14 567 GHz (How et al. 1992), the attenuation coefficient is found 
to be about 0-9dBm _1 for R o %30mf3. The metal loss is the same as for 
conventional waveguides. 



(a) (M 

Figure 5. Distribution of surface current for the TE 10 mode on the walls (a) y = 0, ( b ) y = b. 
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waveguide^ 3 WGCSW,° waveguide 'SkSS ™' l ? = T ^ CWG ’ ^ cntio nal 
walls. 8 6 W,th two conventional and two superconducting 


2.4. Phase velocity, group velocity and phase dispersion 

her^range^tvhlc^detemtne^fiel^Mn^^T' 11 ^ 6 ” 1 penetration dc P th in ‘he tera- 
frequency-dependent skin depth as for normal 0 " ^ thC *2?*®™*’ rather than a 

conductors introduce practically no dispersion into a" 010 ”' ^ mCanS that super * 
in the WGCSW, if any is due Jthe r, ™ . a "“crowave circuit. Dispersion 

by the two broad walls Dispersion cause^h^h ependence the attenuation caused 
velocity is almost the same^ 

and the group velocity * of the TE‘° * 



(24) 


(25) 


wnere c 


is me speed of light in free space. 

dom^na'nTTE^modet smalle^tereatert'lb “ P ?!,?' * y) ° f the WGCSW for the 
waveguides with the same dimensions f * ° Ut than thal of conventional 
for conventional ^““ y is ™*r <l»» .ha, 

geometry is small In summ.rv £. v ld ' d,! P' rs '°" Produced for the same 

WGCSW is reduced drama, ica,7 y ' "» 
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2.5. Maximum transmitted power (dominant mode) 

The maximum transmitted power in the WGCSW is given by 

p “-£ |£ - i ['— (t)T (26) 

where rj 0 is the intrinsic impedance of the medium inside the waveguide, and the 
mangitude of the electric field intensity |£ m „| is given from (15c) as 

l £ »»«l = |tfol^jp (27) 

In the Meissner phase, the maximum electric field intensity and power in the 
WGCSW are smaller than that of conventional waveguides (29kVcm~ l and 1 1 kW) 
because H 0 is very small. For example, they are about 60 V cm -1 and 26 mW for 
H ci % 77mA cm’ 1 , and 3kVcm~\ 64 W for H Ql ^3*85 A cm -1 , respectively. It is 
clear from (26) and (27) that the transmitted power and maximum electric field are 
dependent on the London penetration depth (k c cc\//. i ). The larger the London 
penetration depth, the larger is the transmitted power and maximum electric field. It 
should be noted that the penetration depth also increases with temperature, which 
means that more power is carried by the WGCSW. 

3. Conclusions 

A theoretical analysis for a waveguide with two conventional and two supercon- 
ducting walls has been presented. It is based on high-T c superconducting electro- 
magnetic theory and Meissner boundary conditions. The important implications of 
this paper are: 

(a) the bandwidth of the WGCSW is greater than that of conventional metallic 
waveguides with the same dimensions; 

(b) the tangential electric field and normal magnetic field of the TE 10 mode are 
finite on the inner surfaces of the walls (.x = 0 and x = a), which is in 
agreement with Wang et al. (1994); 

(c) the magnetic field lines are not continuous in the WGCSW, which is the 
result of Meissner effect; 

(d) the surface supercurrents are very large and the surface normal currents are 
very small, so the attenuation coefficient of the WGCSW is smaller than that 
of conventional waveguides by about 10 2 to 10 3 times for different surface 
resistance and penetration depth — the main loss mechanism is metal loss due 
to the two broad walls; 

(e) in the Meissner phase, the maximum electric field intensity and transmitted 
power are smaller than those of conventional waveguides because the critical 
field H cl is very small; 

if) dispersion is small compared with conventional waveguides; 

(g) the wave impedance Z h of the WGCSW has the same expression as for 
conventional waveguides. 
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Scattering by the transition junction between a conventional and a 
high-temperature superconducting waveguide 

RAJ YALAM ANCHILI+. ZHENG AN QIU+ 
and YLN-CHL WANG + 


TE 10 mode scattering during the transition Irom a ^conventional to a >uporeondu^ 
me' rectangular waveguide is investigated based on the Meissner model Thv 
reflection coefficient, input impedance and equivalent inductance are calculated i 
cm-uave bands when the TE I0 mode is incident on the discontinuity, from the 

— A . • r\ «■* \ 11 1 \ llH/> 




1. Introduction 

Recentlv there has been considerable progress in the development ot supereon- 
d uc 1 i n a devices in the microwave and millimetre-wave band, tor example, low - an 
high- temperature superconducting microwave filters, resonators antennas, phase 
Afters etc. The transmission properties of high-T c superconductor wuvegu dc 
(HTSWG) have been studied recently (Wang ei al. 1994). In that paper the 
propagating TE modes in HTSWG are obtained based on Meissner boundary 
condition (i.e. cH.jcn- H.J>. L = 0 on the walls, where H : is the axial component o 
magnetic field. A L is the London penetration depth). . . . 

Since a HTSWG must be connected to its external * /^TSWG 

investigate the EM transition from a conventional waveguide (CWG) to HTS\ G^ 
The physical configuration of the problem under study is shown in Fig. 
turned that the incident wave is the dominant mode, namely. TE l0 mode. 


Y 



Figure 1. Wave scattenng from a rectangular CWG and HTSWG junction at -- = 0 
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travelling from the CWG (z <0) to its HTSWG transition interface (at r = 0), in the 
positive r-direction with exp(j ajt) time dependence. Although the dimensions of the 
CWG and HTSWG are the same, the higher-orders mode TE Pq (p> l y q> 1) (Wang 
et al. 1994) are excited near the transition interface in the HTSWG section (r>0) 
because of the different waveguides. The CWG dominant mode TE l0 is reflected and 
the higher modes TE mn (m#l,n^0) are excited in the CWG as a result of the 
different characteristic impedances of the dominant modes in the HTSWG and 
CWG. Analysis of HTSWG is given briefly in §2. In §3, we consider the reflection 
coefficient, characteristic impedance and equivalent inductance due to the HTSWG 
transition when one, two. three or four modes exist in the HTSWG. 

2. Analysis of high-T c SC waveguides 

A theoretical analysis for HTSWG has been presented by Wang et al. (1994) 
in details based on HTSC electromagnetic theory and the Meissner boundary 
condition. It can be shown that the field component expressions for the dominant 


mode TE 11 in HTSWG are given as follows: 

H : = H 0 (B 0 cos k x x + sin k x x){C 0 cos £ y v + sin k y v) (1 ) 

H x =^X* H 0 (B 0 sin k x x~ cos cos k y y + sin k y y) (2) 

K 

Hy^^rr H o< B o cos k^x-ysin ^x- T )(C 0 sin k y v-c os k y y) (3) 

kc 

i aju n k v 

E x = — — — - H 0 (B 0 cos k x x + sin k x x)(C 0 sin k y v — cos k y y) (4) 

k c 

^r J "TT i H o(B 0 sin M + cos /c ;r x)(C 0 cos & >v v + sin k y y) (5) 

k c 


where fi 0 is the permeability of the material which is assumed to be that of free 
space; H 0 is an unknown coefficient depending on excitations; B 0 = k x k L and 
C 0 = k y /. L are proportionality parameters, but they are different for different modes 
and penetration depth. The wave factor in the form of exp(jcu/ — yz) is assumed; k c is 
the cut-off wavenumber, the real part a of the propagation constant y is the 
attenuation constant, the imaginary part P = 2n!k g is the phase constant, k g is the 
waveguide wave length and cu is the angular frequency; k x and k y are the 
wavenumbers of the TE 11 mode along the x and y axes, respectively. Note that the 
superscripts 1 1 for the dominant mode in the HTSWG have been suppressed in these 
and following expressions for convenience, unless stated otherwise. It is worth 
noting that when considering the HTSWG, because the non-integer mode indexes 
are determined by solving two transcendental equations of dispersion relations 
(Wang et al. 1994) which are derived from the Meissner boundary condition on the 
four walls, in order to avoid any risk of confusion with the commonly used TE W „, we 
use superscripts to designate the modes in HTSWG, i.e. TE M . Here the pair pq bears 



[i \n\'‘juiclc junction s wittermu 


no 


relationship to the wavenumbers. We choose p-\.q 


— ) for the dominant mode 


only for convenience. 

The wave impedance Z M of 

£ 
H 


ihe TE™ mode in the HTSWG is defined as 


E v t 

Z na =.;=- H - jj« 


< > 


whort 




3. Reflection coefficient and characteristic impedance . . = l) 

For convenience. -de 

f° n m'^T?' cul-ofTby «*.* chooying some opera,, ns 
r n " The cut-off modes required for matching the boundary conditions at - • 

frequency. The cut on mo m . TE wa ve ot unit amplitude 

are TE m „ modes in the . r and all the hiaher-order modes are set 

considered. Its reflected wave as amp 1 unknown amplitude a' ‘ 

.« astlwi ,he umgenhal elecric field £. and 
magnetic field H x must be both continuous, i.e. 

x x ' rriTT M7Z nil 

(l+D -sin-.v+ T 1 sin--.* COS jy 

v 71 a m - 1 n = 0 aK cmn 

= a * ( - B 0 sin Arr * + «* cos *>" - v + sin *" • v) <8 ' 

(krv- 


fi l0 a . * f f 

(r _ l) _ r s.n-x+ I I ^ 


mir 

. sin — x cos — v 
a b 


sin Efx-cos A^KCS" cos ^V + sm At^v) <*> 

" (kn l 

where r means ,ha, when r»-l and I »-0. the — "V ^“an” 

“ b0 ' h s,d ” 01 

each of the two equations we obtain from eqns ( 


h n+r) + h{nm)a m „= f a M F(nm, pq) 

° p = i 


( 10 ) 


. . \ x- x /7W 

M£zl> + ^1 I s Hmh.OT) 


[11] 


P = 1 4 “ 1 
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For the HTSWG, when using a similar orthogonality property (see the Appendix 
for the proof), one can obtain from eqns (8) and (9) by multiplying the HTSWG 
wave functions and integrating on both sides of each of the two equations 

g 0 (pq){\ + n+ y f a mn g(rnn, pq) = a pq fipq) (12) 

m = 1 n = 0 


r - 1 ) ^ j, j: 


g(mn. pq) = - ~fipq) 


• 10 


m — 1 n = 0 ~mn 


(13) 


where r 10 is the wave impedance of dominant mode TE 10 . z m „ and a mn are the wave 
impedances and amplitudes of higher-order modes in the CWG; z pq and a pq are the 
wave impedance and amplitude of mode TE pq in the HTSWG, respectively. Other 
quantities are known and given below- 


h 


o — 


a b 


(14) 


h(mn) = 


nhm 

4 HI 


(15) 


F(mtu pq) = 


mnk pq 


+ (- l ) m (B p 0 q sin k pq a - cos k pq a) 


-(W 


1 + (- WCS" sin k”b - cos k pq b) 


, [nn 

(W-( y 


(16) 


a 


1 + cos k pq a — fig* sin k pq a 


- ~{k pq ) 1 2 

a i 


1 -cos k pq b + C pq sin k"b~ 




(17) 


mnk pq 
g(nm , pq) — 


1 + cos k pq a — B% q sin k™a 


1 — cos kfb + C g’sin k pq b 

ZiZZfj 


(Ik) 


f(pq)= 


abk? 

(4 krT 2 


(19) 


In order to solve ( 1 0)-( 1 3) for input impedance and reflection coefficient, etc., we 
use the iteration method according to the order of one, two, three or four modes in 



the HTSWG. Here, only the expressions of one mode are given at the transition 
interface r = 0 (see Fig. 1). The input impedance is 


and the reflection coefficient is 


-1 1 



where Z 11 is the wave impedance of the dominant mode TE 11 in the HTSWG. The 
magnitude and argument of the complex reflection coefficient T. with incidence from 
the'cWG. are plotted in Fig. 2 as functions of frequency when one. two. three or 
lour modes are considered. All magnitude and argument curves show similar 
behaviour. Note that the reflection coefficient magnitude T approaches unity when 
the frequency approaches the cut-off frequency. Above the cut-olT frequency a 
decreasing value of T implies a real power flow into the HTSWG As more modes in 
the HTSWG are considered. |T, and LY converge to the dashed curves. Since the 
characteristic impedances in the two different waveguides are quite different, this 



1.94 1.96 1.98 2.00 


f/tu 


Figure 2. Magnitude and argument of reflection coefficient T as a fu " c ' ion ^ ^ u ^ ncy 
(a = 2b = 2-29cm. /. l = 3-6x 10" ' m. / l0 = 6-55 GHz. / = 126GHz. BW-6 fiSGHz.) 
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type of variation of T is similar to that described by Safaui-Naina and Macphire 
(1982). In Fig. 3, Z m is plotted against the frequency. It should be noted that Z in and 
T are both real quantities when only one mode is considered in the HTSWG 
Similarly, the dashed curves are the results when sufficient modes in the HTSWG are 
considered. It is estimated that 20 modes in the HTSWG will give exact results. 

As shown in Fig. 4. the equivalent inductance decreases with frequency. This is 
due to the superconducting nature of the HTSWG and is independent of the 
dimension-variation type of discontinuities. 

In conclusion, it is found that the transition from a conventional to supercon- 
ducting waveguide gives rise to large reflection { ^09 near j\ {) and decreasing almost 
linearly with frequency). It is suggested that transformers or other devices for 
matching are required. 


Appendix 

Proof of orthogonality of the HTSWG wave functions 

Let u and r (u¥^v) be arbitrary eigenfunctions in the HTSWG obeying the 
Helmholtz equation and the Meissner boundary conditions. We obtain 



Figure 3. Input impedance Z, n = R in +jZ in as function of frequency, (a — 2b = 2*29 cm, 
^ = 3-6 x lO“ 7 m ./ 10 = 6-5GHz,/"= 126GHz, BW = 6 05GHz.) 
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L(mH) 

t 

6.0 -! 
5.0 — j 


4.0 - 



• two modes 
■ three inodes 
▼ four modes 


f /f ;n 


1.9: 


1 . 94 


1 .36 


. 3 8 


Figure 4. Input inductance as function of frequency u - Zb - 2^ '•'m. ' >° m 

h f.„ = 6A5GHz. /"'= 12 6GHz. BU = 6 OiGHz.) 


V J-w + k~u = 0 


i_ w = 0 at .v = 0. x=a: v = 0. y = b 
cn /. L 


V t i' + k 2 t* — 0 


— — J- i? = 0 at x = 0, a = a ; y = 0, y = 6 
cn a l 


then 


l / CW Ct - 5 » . ■> . 7 . 

<b l*s _ “sl +( ‘‘ 


ur d.v dv = 0 


By using Green's theorem, it follows that 

b [t>(V \ + kf)u - u ( V f + k; )r] d.r dy = 0 

From the Meissner boundary condition, we obtain the following relation. 


f: 


(At) 
(A 2) 
(A3) 
(A 4) 

(A 5) 

(A 6) 


<t I' | mi' d.v dv = 0 (AT) 


-j rb 

o Jo 


that is 




(A 8) 


(A 9) 
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Meissner model of superconducting rectangular waveguides 

YEN-CHU WANG+, ZHENG AN QIU+ and 
RAJ YALAMANCHILIt 

The propagation properties of the TE pq mode and their dispersion relations tn a 
high temperature superconductor (HTS) rectangular waveguide using the Meissner 
boundary conditions on all four superconducting walls are presented. In addition 
to recovering some previously known results, w r e have obtained some unique and 
novel results: the dominant mode for HTSWG is the TE 11 mode (instead of the 
conventional TE 10 mode) and the tangential electric field and normal magnetic 
field for the dominant mode TE L1 exist on the surfaces. The expressions for 
electromagnetic components, surface current, attenuation coefficient, the maxi- 
mum transmitted power, dispersion and wave impedance are also presented. The 
attenuation is found to be much smaller than that of the conventional waveguides 
and the dispersion is negligible. The Meissner model presented here has been 
proved to be valid and powerful for analysis of superconducting waveguides. For 
the first time, we have shown that the Meissner boundary condition is a key 
element in boundary value problems for superconductor electromagnetics. 


1. Introduction 

Since the advent of high temperature superconductors (HTS), many supercon- 
ducting electronic devices have been sought at liquid nitrogen temperatures (Nis- 
enofT 1985, Van Duzer and Turner 1981). Incorporating HTS materials in guided 
wave systems at microwave and millimetre wavebands is very promising for high 
speed digital (Hilbert et al. 1989), high frequency analogue (Hammond et at . 1990), 
waveguide and transmission line applications in microwave and millimetre-wave 
devices and circuits (Richard et al. 1992, Winters and Rose 1991, Pond and Krowne 
1988). Some novel applications of HTSC transmission lines and devices have been 
given by Yalamanchili et al. (1992) and Heinen et at. (1991). 

The advantages of using HTS at high frequencies include: firstly, very small 
losses, which means low-attenuation and tow-noise: secondly, very small dispersion 
up to frequencies of several tens of GHz (Kown et at. 1987); thirdly, smaller device 
dimension; and finally the propagation delay time can be greatly reduced because of 
the smaller size and the shorter interconnects (Ei-Ghazaty et al. 1992). 

Experimental and theoretical results have shown that HTS waveguides 
(HTSWG) and transmission lines exhibit significantly better performance than their 
identical metallic counterparts. This is because metal transmission systems have 
ohmic losses at microwave and millimetre-wave bands at which the conventional 
waveguides (CWG) have extremely high signal attenuation (l0 4 dBkm~ l at 
200 GHz (Microwave System Design's Handbook 1987)), and therefore are not 
practical for transmission except over very short distances. The attenuation problem 
can be virtually eliminated through the use of HTSWG and superconducting 
transmission lines. 
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The superconducting waveguides have been studied previously for low T c 
superconductors I < IS k)l Alaux and Wybouw 1976. Rohner 1978. Zepp et al. 1977, 
Fabre et al. 1981). The modelling of the propagation process in HTS transmission 
lines has been dealt with by many authors, using both the two-fluid model and the 
Mattis- Bardeen theory (Dykaareru/. 1988, Lee et al. 1988. Lam et al. 1992. Lee and 
Barf knecht 1992). Hiah-T, (T,>30K) superconductors have different properties 
from low-T superconductors. Vhis paper describes the propagation properties of 
transverse electric (TE) modes in HTSWG based on HTS electromagnetic theory 
and the Meissner boundary condition consistent with the two-fluid model. We will 
describe a direct approach for obtaining the general electromagnetic held component 
expressions, surface current distribution, maximum transmitted power, critical 
current and electrical held for breakdown, attenuation coefficient, dispersion proper- 
ties. and wave impedance and compare them with the CWG. wherever possible. Our 
research conhrms that the dominant mode in HTSWG is the TE' 1 (see below) mode 
instead of the TE 10 mode. Another new property of HTSWG is that the tangential 
electric held and normal magnetic held of the TE“ mode are non-zero on the 
surfaces of walls This property is quite different from the idealized CWG in which 
the tangential E and normal H on the surfaces are equal to zero. 

So far the applied physics and engineering researchers' main interest in the high 
frequency electromagnetic properties of high T. superconductors has been the very 
low attenuation of the EM waves and they have rarely been interested in other, 
usually more intrinsic. EM properties. This research opens up a new chapter by 
emphasizing one of the most important EM properties ot the superconductors, i.e. 
the Meissner effect. It can be shown that the well-known boundary conditions for 
good or perfect conductors are no longer valid and they have to be modified to 
account for the Meissner effect or the diamagnetic property. The new boundary 
condition (called the Meissner boundary condition) gives rise to different field 
solutions hence new waveguide characteristics such as bandwidth, attenuation, 
current distributions, and others. This model is not only applicable to the waveguide 
but also to all microwave and millimetre wave engineering problems as long as 
superconductors, instead of normal metals, are used. 

In §2. the most important dispersion relation of HTSWG is derived from the 
wave equation and the Meissner boundary condition. Expressions for electromag- 
netic field components and surface currents including the surface superconducting 
current J„ and the surface normal current J,„ based on two-fluid theory are given. 
Various parameters of HTSWG are analysed and compared with CWG. The 
conclusions are given in §3. Finally detailed derivations of some formulae are given 
in the Appendices. 


2. Theory 

2.1. Wave equation and numerical solutions oj cut-off wavelength 

In this section, we analyse a HTS rectangular waveguide with all four walls 
superconducting. A rigorous formulation based on Maxwell’s equations and the 
Meissner boundary condition is used to obtain the unique electromagnetic field 
properties in the HTS superconducting waveguide. 

The governing wave equation (Helmholtz equation) and the Meissner boundary 
condition for the longitudinal magnetic field H. of the TE mode in an air-filled 
HTSWG are as follows: 
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V\H : + k'-H z = 0 (') 


l JL-Lh. = 0 at v = 0. x = a: y=0.y = b (-> 

cn r. L 

where the operator Vf is the Laplacian operator in the transverse plane (i.e. vv-plane). 
Propagation in the r direction is assumed. The parameters a and h are the HTSWG 
width and height, respectively. k c is the cut-otT wavenumber. k L is the London 
penetration depth, tor HTS, its values vary from ICT 8 m to 10'" m (Burns 1992). This 
parameter, which is a measure of the distance of magnetic field penetration into the 
superconductor, is very important not only because it can provide information about 
the fundamental mechanisms for superconductivity in materials but also because ot its 
sensitivitv to the quality of the superconductor near its surface. We call (2) the 
Meissner boundary condition. It accounts for the important Meissner elTect in which 
the fields do not vanish abruptly from the surface of a bulk superconductor, rather, 
they decay exponentially into the bulk. The penetration depth r L is the characteristic 
decay length of the magnetic field into a superconductor. It is obvious that when 
y . the second term of the Meissner boundary condition approaches zero, the 

above equations reduce to the boundary condition (cH.,cn)=0 for the CWG. 
Therefore, the first term of (2) is related to the properties of the idealized CWG and the 
secondary term is related to the properties of the HTSWG. Thus the Meissner 
boundary condition can be taken to be of relevance to the two-fluid SC model. 

We assume that the general solution of the governing equation is as follows: 


H. = (A 0 sin k x x + A l sin M)M 2 cos k,y + A 3 sin k y y) (3) 

where 4 0 , .4,. .4,. A } are the unknown coefficients and none of them are equal to 
zero. The k x and k y are the unknown wavenumbers along the x and y axes, respectively. 

Substituting (3) into the Meissner boundary conditions, we obtain two sets of 
equations for k x and k y 




k x sin k x a — - 


cos k x a\ A ( , , sin k x a\ 


(4 b) 


1 


Ai-k,A 3=0 


(5 6T) 


cos k z a\ , (. . ,sink x a\ 

k y sin k y b t- 5 -) A 1 -lk,cosk y b + — -r— J A 3 -0 


(5 b) 


For non-trivial solution condition, the two determinants of coefficient matrix 
have to be equal to zero 
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-it. 


= 0 


cos k x a 

k \ sin k.u — — - — 
'■l 


sin k x a 

— k x cos k x a — 


( 6 ( 2 ) 


-A\. 


cos A v /> 

A\. sin kh 

/-i 


sin k h 
' A\, COS A’ y 6 : 


= 0 


(6 h) 


Therefore, the two important dispersion relations in two transcendental equa- 
tions are obtained from the above determinants as 


tan k x a = 






(7) 


tan k y b = 


^ L k y 
(Ml) 2 - 1 


( 8 ) 


The roots of these two equations describe the dispersion relations as a function of 
the London penetration depth and dimensions for HTSWG. A graphical method is 
used for solution as shown in Fig. 1 from which the cut-off features will be obtained. 

Substituting (4u) and (5 li) into (3). we have (see Appendix A) 

H. = A 0 A, | cos k r x MM Sm kxX ) ( C ° S ky} + k~T Sm ky} ') (9U) 

or 

H z = A l A i (k x /. L cos k x .x + s\n k x .x)(k y /. L cos M + sin k y v) (9b) 


It is obvious that k x ¥=0 and k y # 0 from (9 h), otherwise the solution becomes trivial. 
Our dominant mode is designated to be TE 1 1 . Thus the dominant mode is no longer the 
TE 10 mode. This important conclusion is completely different from that of CWG in 
which the TE 10 mode is dominant. The above formula, however, does cover that of 
CWG in which H. = H 0 cos k x .x cos k y y when k L — x as evident from (9 u). 
Substituting (3) into (1). we have 


therefore, the cut-off wavelength /. c is obtained as lollows 




2 n 

k 


(in 
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Modes 

d=lmm 

a = 2 mm 

a — 4 mm 

a = 8 mm 

TE" 

0-895 (2 0)+ 

1-790 (4 0) 

3-579 (8 0) 

7158 (16 0) 

TE :i 

0-08 (10) 

1415 (2 0) 

2-829 (4 0) 

5 659 (S O) 

TE n 

0-555 (0-804) 

111 ( 1 TS^> 

: :i9 (3-5'si 

4 438 (7 155) 

TE !: 

0-486 M)-7()7) 

0-971 (1414) 

t -941 (2-828) 

3-882 (5-65") 

TE :: 

0-448 (.0 485) 

0 895 (0 9"0) 

1-790 ( 1 940) 

3-579 (3-881) 

TE' J 

0-520 1 0-447 ) 

0-658 1 0-894i 

1316 i l-~S9) 

2 631 (3 5"8i 

+ The lie a re 

s ,n parentheses are the corresponding cut-off wavelengths tor the 

first six TE„„ modes in 

CWG 





Cut-off wavel 

enmhs for some TE P< * modes u L = 10 

m. a = 2b: for /. L = 

10“ 3 m. there is little 


difference between 

the values). 



2.2, Field components and surface current distributions 

The field component expressions for the TE mode can be determined Irom the 
axial magnetic field H z by means of the relations between the transverse and axial 
components as follows. 

From (9 b) % let B 0 = /. L k x , C 0 =/. L k y and = for simplicity. Then the 

magnetic field axial component H : becomes (see Appendix A) 

H Z = H 0 (B 0 cos A* x .v-i-sin k x .\)(C 0 cos Ar v v-r-sin k y y) (12) 

where H 0 is an unknown coefficient depending on excitations. B 0 and C 0 are 
proportionality parameters, but they are different for different modes and penetra- 
tion depth. The wave factor in the form of exp (jtur-yr) is assumed. The real part x 
of the propagation constant y is the attenuation constant, the imaginary part /? is the 
phase constant and tu is the angular frequency. 

Other field components can be expressed as 

i 0 k 

= H 0 (B 0 sin k r x - cos k x .x)(C 0 cos k y v + sin k,y) 

k; 

H y = ^ H 0 (B 0 cos kyx + sin k y .x)(C 0 sin k y y- cos k y y) 

E x H 0 (B 0 cos k x .x + sin k x .x)(C 0 sin k y y - cos k y y) 

k c 

H 0 { — B 0 sin k x .x+ cos k x .x){C 0 cos A-V + sin k y v 

k“ 

where p 0 is the permeability of the material, assumed to be that of free space. 

The wave impedance ZP of HTSWG is defined as 

^ Hy Hy r 


(13 a) 
(13 b) 
(13 c) 
) (13 d) 


( 14 ) 





1 156 


Yen-i.hu ll'uny cl al. 


\1 odes 

a = 1 mm 

u — 2 mm 

c/ = 4 mm 

a =* 8 mm 

TE" 

0895 (2 0)+ 

[ -790 (4 0) 

3 579 (8 0) 

7 158 (16 0) 

TE- 1 

0-08 (10) 

1415 (20) 

2 829 (40) 

5 659 (8 0) 

TE 11 

0 555 (0 894) 

111 ( l '89l 

2 219 (3-578) 

4 438 (7 155) 

TE ,: 

0 486 (0 707) 

097| (| 414) 

1941 (2 828) 

3-882 <5 657) 

TE : ’ 

0 448 (0-485) 

0-895 (0 970) 

[ 790 ( l 940) 

3-579 (3-881 ) 

TE ,J 

0 529 (0 447) 

0-658 (0-894) 

1 316 (1-789) 

2 631 (3 578) 


+ The riiiures in parentheses are the corresponding cut-otT wa\e!cngths for the rirst six TE^ modes in 
CWG 


Cui-olT w avelengihs for some TE 79 modes {/. L = 10 m. a = 2b: tor /. L — 10 * m, there is little 

difference between the values). 


2.2. Field components and surface current distributions 

The field component expressions for the TE mode can be determined from the 
axial magnetic field H. by means of the relations between the transverse and axial 
components as follows. 

From (9 b), let B 0 = /. L k x , C 0 = k L k y and H 0 ~A l A 3 for simplicity. Then the 
magnetic field axial component H. becomes (see Appendix A) 

H. - H 0 (B () cos k x x 4- sin k x x)(C n cos li,y* sin k y y) (12) 

where H 0 is an unknown coefficient depending on excitations. B 0 and C 0 are 
proportionality parameters, but they are different lor different modes and penetra- 
tion depth. The wave factor in the form of exp (jcur — yr| is assumed. The real part a 
of the propagation constant y is the attenuation constant, the imaginary part /? is the 
phase constant and a> is the angular frequency. 

Other field components can be expressed as 

H x - 3 'TT w o< s o sin k r x- cos k z x)(C 0 cos k v y + sin k y y) 

H= 3 -r-r H 0 (B 0 cos fc^x + sin k x x)(C 0 sin k i -cos k y y) 

K ; 

E x ~ — Tf- 1 W 0 (S 0 cos kfX + sin MKQ sin k r y- cos k t y) 

E y h 0 ( - B 0 sin k x x + cos k^xHCo cos k y y + sin k y y ) 

where is the permeability of the material, assumed to be that of free space. 

The wave impedance of HTSWG is defined as 

h H y h x r 


(13a) 
(13 6) 
03 c) 
(13 d) 


( 14 ) 
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Figure Distribution of electric field E x and E y for TE 11 mode in a cross-section, (n) E x . 
* (h) E y . 


The instantaneous field expressions for the dominant mode TE 1 1 in HTSW G can 
be obtained as 

H. = H 0 (B 0 cos k x x + sin k x .x) ( C 0 cos k y y + sin k y y) cos (tot - pz) 0 Su) 

H = ft HJ B 0 sin k x x - cos k x)( C 0 cos k y y + sin ky) sin |wf - pz) ( ' 5 b) 

H v =- ft H 0 ( B 0 cos k x x - sin k x x) ( C 0 sin ky - cos ky) sin M -pz) (15 c) 

E = H 0 (B 0 cos M + sin k x x)(C 0 sin ky -cos ky) sin {au-pz) (\5J) 

x k ; 

E = H 0 ( - fl 0 sin (fc x .v + cos jfc,x>( C 0 cos + sin ky) sin (cur -pz) (15c) 

> k: 


where k P , and A, are the cut-off wavenumber, the phase constant, and the waveguide 
wavelength of the TE 1 1 mode, respectively. k x and *, are the wavenumbers of the TE 
mode along the * and y axes, respectively. Note that the superscripts 1 1 for the 
dominant mode have been suppressed here and below unless otherwise stated. 

The above equations (15 </) and (15 e) demonstrate obviously that due to B 0 and 
C 0 being very small, the space distribution of the electric field £ y of the TE mode is 
proportional significantly to cos kyc sin ky, which means that the tangential electric 
fields exist on the surfaces (x = 0 and .x = a) of HTSWG and has maximum at v *b - 
The space distribution of the electric field E z is proportional mainly to 
sin M cos /c y, which also means that the tangential electric fields exist on the 
surfaces (v = 0 and y = b) of HTSWG and are a maximum at xvai 2. The unique and 
novel property of the finite electric field tangential components on the surfaces of 
HTSWG was obtained first, which is completely different from that of a perfectly 
conducting CWG since, for H t * 0 and for a finite surface impedance at microwave 
and millimetre wavelength, the tangential components of £ on the walls are finite. 
For details of the field distributions see Fig. 2. 
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y 




y 



U*) 

Figure 3. Distribution of magnetic field H x , H y and H. for TE U mode in a cross-section. 

(a) H x : {b ) H y : ic) H z . 


The above equations (15 a), ( 15 /?) and (15 c) demonstrate obviously that as B 0 
and C 0 are very small, the space distribution of the magnetic field H y of TE U mode 
is significantly proportional to sin k^x cos k y y. which means that the tangential 
magnetic field is very small on the surfaces (x = 0 and x = a) of HTSWG and normal 
magnetic fields exist at surfaces of v = 0 and y = b. The space distribution of the 
magnetic field H x is proportional mainly to cos/r^vsin/y, which means that the 
tangential magnetic field is also very small on the surfaces of y = 0 and v = b of 
HTSWG and the normal magnetic field exists on the surfaces of x = 0 and x = a. On 
the other hand, due to H : x sink r x sink the magnitudes of the tangential 
magnetic fields are negligible on all four surfaces. Details of field distributions are 
given in Fig. 3. The above results demonstrate that the magnetic field lines are not 
continuous in HTSWG (Orlando and Delin 1991), which is in agreement with the 
Meissner effect, but the magnetic flux density B is continuous in and out of HTSWG 
walls because B = /i 0 (H +M) (M is the magnetization density in the walls and is not 
equal to zero). This also is a unique property in HTSWG, which does not exist in 
CWG. 

The London equations and the two-fluid model are used to investigate the 
current distribution. The total conduction current density consists of the supercon- 
ducting current density J, and the normal electron density J„ based on the two-fluid 
model. This model assumes that the conducting electrons in the HTSWG walls are 
divided into two categories. We use a different method to describe their distribution. 
The surface superconducting current due to superelectrons is calculated based on 
the London relation and appropriate boundary condition first. Then the normal 
surface current due to normal electrons and the losses arc calculated using the 
approximate resistive boundary condition and the perturbation approach. 
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The superconducting current J, in four isolated walls model in which each wall is 
considered as a simply-connected superconductor using thin film substrate and can 
be given as 

j 5= -i A < 17 "» 

This equation is called the London relation, where A = «„/.[. A. in units ot 

Tesla-meter, is the vector potential in HTSWG. Let the vector potential inside the 
air-fiiled hollow region in the waveguide be Aj It satisfies the following equations. 

V x Aj = B H 7 ^> 


Thus London relation guarantees that the vector potential is a Teal field) that is 
defimtelv specified. Equation (17c) is called the Lorentz gauge, which gives rise to 
the relation between electric field E and vector potential A a in the air region. 

Before solving the current distribution problem, it is important to determine the 
boundary conditions at a boundary between the two media (here the superconduct- 
ing medium and the air). The boundary condition on the tangential components ot 
fields (J s and A a ) for the simply-connected SC is 

I = - ‘ A (IS) 

J \t ^ 


where the subscript 'at' denotes the component tangential to the boundary in the air 
or the applied field region, the subscript 'st' denotes the component of supercurrent 
density tangential to the boundary in the HTS. Note that )„ is the bulk current. 
Once A„ is known, the current density K J( . and surface supercurrent density J ss , are 
readily calculated. The solution satisfying (17 b) and (17c) is obtained from B as 
follows 

A a: = ,4 0 (B 0 sin k x .x + COS k x .x)(C 0 sin k y y + cos k y v) < l9tJ) 

A ( S 0 cos k x x - sin k r x) (C„ sin k v y + cos k y y) 

/> 

H 0 (B 0 cos k x .x + sin k x .x)(C 0 sin k y y -cos k y \ ) ( 19 £ 


A ay j} 


jA ° Ky (So sin k r x + cos k x .x)iC 0 cos A. v y-sin k,y) 


h 0 (B 0 sin k x .x - cos k t .x)(C 0 cos L v v + sin k f y) 


k: 


(19c) 


The surface supercurrent distributions on the isolated walls are obtained (see 
Appendix B) as follows 
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y = 0 plane: 

Jtsi — J "I” ^ssr^r 


A 


fifyk. v JA 0 k x 

-pr Hq(B 0 cos & x .x4-sin ^ r x) -f — ( B 0 cos A x .t — sin A*.*) 

A:; p 




— (£ 0 s j n -,-cos Ar^rja. 

A 


(20 12) 


v = b plane: 


J„, = A,A + •/„=*: 


A 


a r» At , . } A r* Ar r 

--^2 tf 0 (£ 0 cos A: x .v + sin Ar x .v) H — — (B 0 cos k x .x-sin k x x) 

A'* P 


4- (£ 0 sin Ar t .v + cos A x .v)a. 

A 


(20 b) 


x = Q plane: 

J J 5 i J ssy^y ^ ss;^: 


A 0 /. l 


PqA j A (yk ~j 

- ^o(C 0 cos A y v + sin k y y) +' — (C 0 cos A y y - sin k y y) ; 

(C 0 sin A: y v4-cos /c y >’)a : 


(20 c) 


.t — a plane: 

Jh( = J ssy^y 2jji-a ; 


A 


H 0 (C 0 cos A y y4-sin k y y) + jAoky (C 0 cos A y y-sin A y v) a y 
A c P J 


.4 A 

4-—^ (C 0 sin A: v 4- cos k v y) a_ 
A 


(20 1/) 


where J„ x , J wy and J„, are surface supercurrent magnitudes on the walls along the 
t, v and z directions, respectively. We consider here that the superconducting current 
decays exponentially in the walls, therefore the current is mainly carried by the 
superconductor surface adjacent to the air in the HTSWG. It is easy to see that the 
surface supercurrent and bulk current are very large because the amplitude of J„, is 
proportional to \ jk L (A = ^ 0 kj) which is about the order of 10 7 . All high-T f materials 
are type II superconductors. In the Meissner phase, we assume that the critical field 
/f cl 3:7-7 to 385 Am* 1 (depending on temperature) for HTS sample 
YBa 2 Cu 3 O y (88 2A, H~c axis) (Wu and Sridhar 1990), therefore H 0 *: 13 to 
660 Am* 1 , then the theoretical values of supercurrent density are estimated to be 
about 10 6 * 8 A cm* 2 which is in excellent agreement with other data (Bums 1992, 
Hemen et al. 1990, Miranda et al. 1991, Levenson et al. 1991). The distribution of 
surface supercurrent of HTSWG is shown in Fig. 4. (It should be noted that the 
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on the x-a plane: 

Jsn = J snz^z ^ ^ sny*y 

H 0 (B 0 cos k K a- sin k x a)(C 0 sin k y y- cos k t y)a. 

— H 0 (B 0 cos A.' t a + sin k z a){C 0 cos A^.i -t-sin k v \ )a, (21 d) 

where J mx . J„ v and J are the surface normal current magnitudes on the walls 
alons the v. v and r directions, respectively. It is easy to see that the distribution is 
the same as that of the surface supercurrents. The losses of HTSWG come from the 
surface normal currrent flowing through the small surface impedance. It is predict- 
able that the losses are very small because factors B 0 and C 0 are very small. The 
distribution of surface normal current of HTSWG is shown in Fig. 4. 

2.3. Attenuation coefficient 

The attenuation (caused by the surface normal current) based on the two-fluid 
model in HTSWG walls can be calculated by the perturbation technique as given by 
(see Appendix C) 

pq Ptj£) o lk; k;P>+ktik; + 2k ;) 

* P m bunjk; 

For the TE 11 mode. a = 2mm, 6= 1mm. k x = Zrt, 4 00 1 7 mm, /c y = 2rr/2 001 mm. 
k c = 2tc 1 -79 mm. substituting into (22). we have approximately 



where/, is the cut-off frequency of the TE 1 1 mode. R, is the surface resistance of the 
HTSWG walls in ohms. The unit of >. L is the millimetre. Note that the attenuation is 
not only proportional to f] because R,x/ 2 (Romanofsky and Bhasin 1991), but is 
also proportional to k 2 L . The relation curve of the attenuation coefficient a versus the 
frequency f is shown in Fig. 5 where R,= 50x 10 “®/ 2 Q and a = 2b = 2 mm. 

The losses in HTS waveguides are generated from two sources: superconductor 
loss and dielectric loss. Typical superconducting waveguides used in microwave and 
millimetre-wave circuits use low-loss dielectrics or free space, so the dominant loss 
mechanism is the superconductor loss. This kind of loss is caused by the current of 
normal electrons in superconducting bulk based on the two-fluid model theory. The 
paired superconducting electrons are dissipationless because they cannot be scat- 
tered without breaking pairs. The microwave surface resistance of superconducting 
films has been found to be as low as the order of micro-ohms. These surface 
resistance values are lower than those for Cu at the same temperature by a factor of 
several tens. Such a surface resistance gives a small power loss. For example, 
R,= 1 16mfl (77it) at 58-6GHz (Heinen et al. 1990). x L =lO- 7 m, the attenuation 
coefficients are found to be about IT x 10 "dBkm l , which is much smaller than 
that of CWG (see Fig. 5). If the penetration depth A L = 3 6x 10' 7 m for the case of 



\fe 


•/ ,/' .\aieauntes 


1 I* 


conduction along the 
14 567GHz (How et 
x 10 -3 dB km“ l . 


c-axis and the surface resistance is equal to 10*3 mH (60k) at 
at 1992), the attenuation coefficient is tound to be 9 


^ 4 Phase velocity, group velocity and phase dispersion 

" Superconductors have a frequency-independent penetration depth to the tera- 
hertz ranee that determines field penetration into the material rather than a 
frequence dependent skin depth as for normal conductors. This means that super- 
conductors introduce practically no dispersion into a microwave circuit. Dispersion 
in HTSWG if anv, is due to the frequency dependence in attenuation which has 
been shown to be'negligible because of the very small attenuation coefficient. The 
dispersion caused by frequency dependence of the phase velocity ,s a most the same 
as for the CWG. The phase velocity i „ and the group velocity r, of the TE mode 
can be shown to be given by 


phase velocity 



I 24 i 


group velocity 



(25) 


where c is the speed of light in free space. 

Our work shows that the phase velocity v, of the HTS waveguide for the 

dominant TE 1 1 is slightly smaller (about 01%) than that of 

euides with the same dimensions and frequency because/, is smaller than that ot he 

CWG. Therefore, the dispersion produced for the same geometry is practically the 



Figure 5. Attenuation due to wall losses versus frequency. %-10' m. 
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same. In summary, the signal distortion due to dispersion m the HTS waveguide is 
reduced dramatically because it has no dispersion due to attenuation. 


2.5. Maximum transmitted power (dominant mode ; 

The maximum transmitted power in HTS waveguide is given by (see Appendix 
D) 


P mil 



(26) 


where r\ Q is the intrinsic impedance of the medium inside the waveguide and £ 
the magnitude of the electric field intensity is given by 


In the Meissner phase, the maximum electric held intensity and power in 
HTSWG are smaller than that of CWG (29kVem and II kW). For example, they 
are about 60Vcm, !3mW for H cl ^77mAcm _l and 3kVcm~ L , 32W for H cl 
^3-85 A cm “ l , respectively. It is clear from eqns { 26) and (27) that the transmitted 
power and maximum E are dependent on the London penetration depth ( k c ~x I /. L ). 
The larger the London penetration depth, the larger is the transmitted power and 
maximum E. It should be noted that the penetration depth also increases with 
temperature, which means that more power is carried by the HTS waveguide. 


3. Conclusions 

A theoretical analysis for the HTS rectangular waveguide is presented. It is based 

on HTS electromagnetic theory and the Meissner boundary condition. The import- 

ant implications of this paper are: firstly, the bandwidth of the HTS rectangular 
waveguide is only 1 27 : 1 which is less than that of a conventional metallic waveguide 
with the same dimensions; secondly, the tangential electric field and normal 
magnetic field of the TE 11 mode are finite on the inner surfaces of walls, which 
demonstrate another quite different property between the HTSWG and the idealized 
CWG; thirdly, the magnetic field lines are not continuous in HTSWG, which is the 
result of the Meissner effect; fourthly, because the surface supercurrents are very 
large and the surface normal currents are very small (the attenuation is very small), 
the attenuation coefficient of the HTSWG is smaller than that of the CWG by about 
lO^-lO 3 * * * 7 times for different surface resistance and penetration depth. Also, the 
dispersion caused by attenuation is negligible. This is also quite different from the 
CWG; fifthly, in the Meissner phase, the maximum electrical field intensity and 
transmitted power are smaller than that of the CWG because the critical field H ci is 
very small; sixthly, the dispersion is negligible compared with the conventional 
waveguide, as there is no dispersion due to attenuation; and finally, the wave 
impedance Z £* has the same expression as the CWG. 

Superconducting waveguides distinguish themselves from normal metal wave- 
guides in the microwave and millimetre-wave bands. Firstly, much lower surface 
resistance can be obtained with superconductors, giving rise to much lower loss. 
Secondly, the superconductors have a frequency-independent penetration depth that 
determines the field penetration into the material rather than the frequency- 
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same In summary, the signal distortion due to dispersion in the HTS waveguide is 
reduced dramatically because it has no dispersion due to attenuation. 


2.5. A fa. xtmum transmit ted power dominant mode 

The maximum transmitted power in HTS waveguide is given by (see Appendix 
Dt 


P max 



(26) 


where is the intrinsic impedance of the medium inside the waveguide and E . 
the magnitude of the electric field intensity is given by 


In the Meissner phase, the maximum electric held intensity and power in 
HTSWG are smaller than that of CWG (29kVcm and 11 k\V). For example, they 
are about 60 V cm, 13 mW for ^ 77 mA cm " 1 and 3kVcm' 1 , 32 W for H c[ 
^3*85Acm~ l , respectively. It is clear from eqns (26) and (27) that the transmitted 
power and maximum E are dependent on the London penetration depth ( k c x 1 /. L ). 
The larger the London penetration depth, the larger is the transmitted power and 
maximum E. It should be noted that the penetration depth also increases with 
temperature, which means that more power is carried by the HTS waveguide. 


3. Conclusions 

A theoretical analysis for the HTS rectangular waveguide is presented. It is based 
on HTS electromagnetic theory and the Meissner boundary condition. The import- 
ant implications of this paper are: firstly, the bandwidth of the HTS rectangular 
waveguide is only 1-27:1 which is less than that of a conventional metallic waveguide 
with the same dimensions; secondly, the tangential electric field and normal 
magnetic field of the TE 11 mode are finite on the inner surfaces of walls, which 
demonstrate another quite different property between the HTSWG and the idealized 
CWG; thirdly, the magnetic field lines are not continuous in HTSWG, which is the 
result of the Meissner effect; fourthly, because the surface supercurrents are very 
large and the surface normal currents are very small (the attenuation is very small), 
the attenuation coefficient of the HTSWG is smaller than that of the CWG by about 
10 3 4 -10~ times for different surface resistance and penetration depth. Also, the 
dispersion caused by attenuation is negligible. This is also quite different from the 
CWG; fifthly, in the Meissner phase, the maximum electrical field intensity and 
transmitted power are smaller than that of the CWG because the critical field H cl is 
very small; sixthly, the dispersion is negligible compared with the conventional 
waveguide, as there is no dispersion due to attenuation; and finally, the wave 
impedance Z% q has the same expression as the CWG. 

Superconducting waveguides distinguish themselves from normal metal wave- 
guides in the microwave and millimetre-wave bands. Firstly, much lower surface 
resistance can be obtained with superconductors, giving rise to much lower loss. 
Secondly, the superconductors have a frequency-independent penetration depth that 
determines the field penetration into the material rather than the trequenev- 
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dependent skin depth characterizing normal conductors. This means that supercon- 
ductors introduce no dispersion into a microwave and millimetre-wave system. 
Based on these unique characteristics, the HTS waveguide is reasonably attractive 
for several applications. 

This paper has reported that the Meissner model presented here has been proved 
to be valid and powerful for a more rigorous analysis of superconducting electro- 
magnetic waveguides than those previously attempted which only considered the 
surface resistance. It should be noted that the Meissner condition has fully 
accounted for the EM property of an SC surface in the microwave and millimetre- 
wave bands although it is only indirectly related to the attenuation or surface 
impedance which is due to the presence of both superconducting and normal 
electrons, as in the two-fluid model. This is the first successful attempt to model the 
propagation characteristics of HTSWG using the Meissner boundary condition. We 
have shown that the Meissner boundary condition is a key element in boundary 
value problems in superconductor electromagnetics without which its intrinsic 
electromagnetic properties will not be complete. 

This paper has discussed only the case of four isolated SC walls (thus four simply 
connected SCs) in order to simplify the calculation of the surface current. In the 
future the case of doubly connected SC waveguide with all four walls as one piece of 
SC will be considered wherein the surface current will be seen to depend on the 
magnetic flux quanta passing through the waveguide. Also the anisotropic resistivity 
effect on the SC waveguide propagation as well as the scattering properties will be 
presented. 

Appendix A 

Axial magnetic field H, 

From (4a) and (5 a), two relations can be obtained as follows 

A 0 = k x A. L A i (A I a) 



Appendix B 

Calculation of surface superconducting current 

y — 0 plane: 

From (18), we have 



1 1 ro 


Yen- chit W ;nj et al. 


J st = J iX =~ A A u 


(B l) 


Therefore, the surface currents can be obtained based on the Meissner effect as 
follows 

:'0 

•/>«=! exp (.» / j Jy 

u - r 

J. * ! *) 

= ' U,A ; H ,(B cos k'.x+ sin k x x) 1 (C, sin Ayr -cos Ayr) exp (y ). L ) d.v 

As- J - r 


. j /. r o 

- 0 ‘ (S„ cos Ay.v-sin Ayr) (C„ sin Ayr + c os Ayr) ex P 0 '/.> dr 

A/i J - / 


x/ t> Ar „ , . 

= - ‘ , cos A: v v ^ sin A: t .v ) 

Ak; 

J A x . . 

{ B ) cos k r x — sin k x x) 

Mi 


'■L 

1 -r ik L k v ) z 

1 - ( /. { k,)~ 


( B 2l 




.^4 


*4 1 u 

-° (S 0 sin Avv + cos Ay.x) I < c o sin Ayr + cos Ayr) exp (,v /. L ) dr 

A J - , 

'l-Co'-'l '% 


-- ( Bq sin + cos k x x) 


1 -H a l A v )~ 


( B 3) 


Because C 0 and A L are very small, we have 


J,„ ' - H ol So cos A,.x + Sin Ay.v) < Bo cos Aye - sin Ayr) ( B 4) 

Ak; A/J 


*4 / 

J s . ^ ( B 0 sin Ay.v + cops A t .v) 

A 


v = b plane: 

-/„*=( hx exp (-<y-A>) /. L )dy 

Js 


(B 5) 


= £ W 0 (S 0 cos A,.v + sin k x x) I (C 0 sin k t y-cos k y v) exp ( -(y-Sj/Aj dr 

Ak: L 


* 'A/A 


(B 0 cos /c x ,x — sin k x x) 


if (C 0 sin k y ) 


-^. 2^1 H 0 (B 0 cos Aye + sin Aye) 
Ak; 


, v + cos k y y) exp ( — (y — b) : k L ) dy 

(B 6) 

(Cq k L ^/-ik y ) sin k v b + (C 0 kik y — k L ) cos k y b 


J ssz = - (B 0 sin k x x + cos k 

A 

.4 


= - Y ( B„ sin Aye + cos k x x) 


l +(/. L k f ) z 

Sy.r) j (C 0 sin Ayr + cos Ayr) exp ( -(y-b),/. L ) dr (B 7) 
1 *> 

( d q/~l — y ) ~b (^O^L^y ^i.) C0S ^ yb 


1 -M/. l A: v ) 


V \2 


>s A: v 6~ 



Meissner model of superconducting rectangular waveguides 


1167 


Because C 0 and X L are very small, and cos k y b ^ — 1, then 

Link. .A r iAnk 

J 


\k* 


JA 0 k x A L 

H 0 (B 0 cos kjX + sin k x x) H ^ — (£ 0 cos k x x — sin AvO 


A/? 


(B8) 


° - (B 0 sin A^x-hcos A r x) 


(B9) 


Similarly, we have J ssy and J a: on the ,x = 0 and x~a planes. 
,x = 0 plane: 


■ r*U'*x u , ; , ■ / \ Bo'-L^x + -‘•L 

•/«„ = TT2 W o (C ° cos V' + sm | + (; t & ~~p~ 


Mp*x 

: Aik; 


I ^ 

- ° - (C 0 cos A: v — sin k 

Ap 


, 'l ~ B 0 /-ik x 

y> 1+U^ *P _ 


(B 10) 


A 

■/,„= (C 0 sin A y y + cos *„>>) 


Because B 0 and are very small, we have 


V *-L+ x ~\ 

L (^l*x) 2 + 1 J 


Ho(C 0 COS A y y + sin (C 0 cos A y y-sin k y y) 


(Bill 


(B 12) 


A X 

J„z « — jf [C 0 sin k y y + cos k y y) 


(B 1 3) 


x = a plane: 

J ssy~ ~ i ^^r Ho(Co COS k y y + sin k t y) + ^ A °^ ~ (C 0 cos A y y-sin k y y) (B 14) 


J 


JSZ 


AqX l 

A/J 


(C 0 sin kyy + cos k^) 


(B 15) 


Appendix C 

Calculation of attenuation coefficient 

Using the commonly accepted formula for power in waveguide electromagnetic 
theory, it is easy to show that transmitted power or the time-average power flowing 
through a cross-section of HTSWG for the TE pq mode may be expressed as: 
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yV'T-i. hu et a 1 . 
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= \ R,.Z? I “ I \h x H:-H,H^ d.vd.v 

J 0 J 0 


'> 2 / . 2 ^ j ** h , 

= xt ‘ H„ : Z; 9 I IB, sin A.\.v-cos k t .x) z {C 0 cos A v .r-sin A\.v)* d.vdv 
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Bo 
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In order to calculate the time-average power loss caused by the surface normal 
current in the walls per unit length, we must consider all four walls. Owing to 
j ( 1 - 4 - 0 )= _j i(i (, =B).J sn (\=0)= - J J v = u), the total power loss is then double the 

sum ot the losses in the walls at v = 0 and v = 0. We have 


P L (-) = -[P L ( 0 + -[Pt( - >lv = 0 


(C 2) 


where 


[P L (r)] y ,o = i j V-O'-OM'd* 




JO 


\H () \ 1 Cq(B 0 sin k x x — cos k^x) 1 
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+ ! H 0 \ 1 C 2 o(B 0 cos k x .x + sin k.x) 2 dx 


x-O- 


= ’ R,\H 0 \ 2 C's 


.j, . 1 zB 2 0 
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Bq i 


• (1 + Bs) 4 - - sin 2k x a+~ sin - k x a 


a , l — Bq . # B 0 . , . 
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l ^ 
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•f sn( X = 0)1 : d.V 




'B 2 k 2 

p- |H 0 | 2 Bo(Q sin £,y-cos/c y j) ; 


+ \H 0 \ 1 B^{C 0 cos fc y j*4-sin A: y > ) 2 


dr 
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- ^ R s \Ho\ m Bo ) jpr 
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4/C>, AC y Jj 


(C 4) 


Owing to stn 2 k x a. sin 2 k x a. sin 2 kjb. sin 1 k y h. B 0 and C 0 being very small, the 
higher small terms are ignored and we have 


ab , /? 2 Zf« aA 
P,=s- IHo! J? 1 rr- = j iHoko 


1 - 


,n g 


(C 5) 


where ZJ* is the wave impedance of TE P * mode. 

Substituting (C3) and (C4) into (C2) and omitting the higher small terms, we 
have 


P l (:)A<*R,(\Ho\C 0 ) 2 \ 1 + 


pl ^ + l -bR,(\H 0 \B 0 ) 2 ^\+^j 


= - rj// 0 I 2 ^£ ^ 


lW + 4{W + «]} (C6) 


Therefore the attentuation constant of TE M mode can be obtained as follows 
(a = 2b): 


, _ P £.(-) _ R . 2 3 kiklp 1 + kt(kl + 2* 2 ) 

2P m to*,**, 2 


(C7) 


For the dominant TE 11 mode, a — 2 mm, 6= 1mm, k x = 271/4 001 7 mm, k y 
= 2rr/2 00 1 mm, k c = 2n/\ -79 mm, substituting into (C 7), we have approximately 


/c 


a n ^2-2x 10 -4 /?^ 


ffM' 
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l - 
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(C 8) 


,/ 


where is the London penetration depth in millimetres. 
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Appendix D 

Calculation of maximum power 



uhere f:„ a N the maximum electric Held tor breakdown. 
From i D 1 ) we have 


tf () 2 = 


, k < 


= E r 


Substitution of i D 2) into iC 5) yields 
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(D 1) 


( D 2) 
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Using the cavity model and Green's functions, equations are derived to express the input 
impedance of a rectangular patch microstrip antenna as a sum over a single integer index. 
It is shown that the computer execution time for calculating the input impedance is greatly 
reduced using the single sum approach rather than the popular double sum approach , without 
sacrificing accuracy. Double sum solutions were found to take at least twice as long to 
perform as those of the single sum solutions. The validity of the theoretical results was 
verified by comparing them with experimental results. 
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Xp x-coordinate of probe location 

y y-coordinate of probe location 
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ko 
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E, 


transverse Laplacian operator 
free space wavenumber 
rectangular-patch (cavity) wavenumber 
cavity modal wavenumber 
(-1 ) m 

z-directed electric field confined within the cavity 
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Notation 

Xp x-coordinate of probe location 

y p y-coordinate of probe location 

V, 2 transverse Laplacian operator 

k 0 free space wavenumber 

k rectangular-patch (cavity) wavenumber 

cavity modal wavenumber 

j (-ir 

E z z-directed electric field confined within the cavity 
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J 

h 


Ho 

Io 

S(x-Xp) 

5(y-y P ) 

^mn 
4^ urn 

a 

b 

e(m) 

e(n) 

Z 

Aw 

Zss 

Zjn,iK 

Z«U.SS 

«x) 

gm(x) 
A,,, B„ 


single-sum cavity electric field for y £ y p 

single-sum cavity electric field for y s y p 

probe (source) current density 

z-directed probe (source) current density 

cavity angular frequency 

modal cavity angular frequency 

permeability of free space 

Amplitude of z-directed current density (source) 

Dirac-delta function used to model x-directed current density contribution at probe location 

Dirac-delta function used to model y-directed current density contribution at probe location 

complex modal constant associated with the z-directed electric field 

modal function associated with the z-directed electric field 

patch length 

patch width 

Neumann number (=1 if m = 0, = 2 if m* 0) 

Neumann number (=1 if n = 0, = 2 if n* 0) 
probe impedance 
double-sum cavity impedance 
single-sum cavity impedance 

double-sum cavity input impedance at probe location 
double-sum cavity input impedance at probe location 
x-dependent single-sum electric field modal function 
x-dependent single-sum electric field modal function 
complex constants associated with ^(x) 
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p real constant associated with fjx) 

H magnetic field intensity 

H, x-component of magnetic field intensity 
f (y) single-sum sinusoidal function for y > y p 

f<(y) single-sum sinusoidal function for y < y p 

C m , D m complex constants associated with f„(y) 

E m , F m complex constants associated with f«(y) 

A,„, A' m complex constants associated with g„,(y) 
y(= Ym) single-sum propagation constant 
a x x-directed unit vector 

H inward-directed unit vector normal to cavity wall 

Q ohms 

1. Introduction 

Several authors have produced fine work concerning microstrip antennas and their various properties, such 
as input impedance, Q-losses, radiation, etc. (Carver and Mink 1981, Howell 1975, Pozar 1982, Agrawal 
and Bailey 1977, Dereyd and Lind 1979, Deshpande and Bailey 1982, Watkins 1973, Lo, et al 1977). 
Many previous works (James and Hall 1981, Lo et al. 1979, and Richards etal 1981) have used the cavity 
model analysis to derive the electric field and input impedance of various microstrip antenna configurations. 
Usually, the plane of the microstrip is viewed as one geometric region, with length a and width b. This 
approach provides quantities which depend on double-summation solutions summed over two integer mdtces. 
Other authors (Chadha, et al. 1981, Alhargan et al 1991 ) have used reduced single summation techniques 
to find Green's functions for microstrip circular disk and annular ring configurations. In this paper, the 
rectangular-patch microstrip antenna is analyzed and theoretical and experimental data are directly compared 



4 

By viewing the plane of the microstrip patch as being two separate regions, I and II, bordered by the 
ordinate of the antenna's thin-probe coordinate, y = y p , the electric field and the input impedance can be 
expressed in terms of solutions which are summed over only one integer index. This single-sum approach 

greatly reduces execution time, as compared to the double sum approach. 

The execution times are compared, and the efficiency of the single-sum approach is demonstrated. 

2. Theory 

2. 1 . The cavity model 

Authors Lo, Solomon, Richards and Harrison (1979, 1981) developed an analysis of microstrip antennas 
by modelling them using a lossy cavity bounded by electric walls on the top and bottom, with magnetic walls 

on the sides. 

Figure 1 is the physical representation of a cavity model. The model applies to thin-substrate microstrip 
antennas whose substrate thickness is s 1% of a free space wavelength. Fields underneath the top conductor 
(patch) and above the ground plane are determined by the cavity model such that the electric fields in the 
antenna structure and the input impedance at the probe position can be determined. 

The cavity model assumes that the magnitudes of the tangential magnetic field components at the edges 
of the patch are negligible, implying a high-impedance condition at the boundary. This condition is equivalent 
to modelling the perimeter of the patch with magnetic walls. The electromagnetic field components within 
the cavity are obtained by solving the appropriate Maxwell's Equations within the cavity, subject to the 
appropriate boundary conditions. It is assumed that the electric fields within the patch are z-directed and 
invariant with respect to z, due to the presence of a thin substrate. 

2.2. Double-sum approach 
2.2.1. Electric field derivation 

The z-directed electric fields under the patch due to a z-directed current source (probe) must obey the 



Electric Well 


5 




inhomogeneous Helmholtz Equation : 
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( V f 2 + k 2 ) E t * ( J • ? ) - o 


( 1 ) 


where, J z is the z-directed current density on the probe. Since such probes are usually very thin with respect 
to the excitation frequency, delta functions have to be used to describe the current density 

J t - 7 0 6(x-xp biy-y r ) (2) 

where, (x„,y p ) is the coordinate of the probe position on the patch. Assuming E z can be expanded in terms 
of a Generalized Fourier Series, 


m-0 n * 0 


(3) 


where, a„„ is a constant , and Y^(x,y) = E 0 Cos(mnx/a) Cos(nny/b), such that E 0 = V 0 /d is constant, 
d = substrate thickness , and the 's represent modal solutions to the homogeneous form of (1). 
Substituting (3) into (1) , 

E E -O*- * (4) 

m-0 n-0 


where, k^ 2 - (mnx/a) 2 + (nnx/b) 2 . Multiplying both sides by T BM *(x,y) and integrating from x - 0 
to x = a, and y = 0 to y = b, over dS = dxdy , the z-directed double-sum form of the electric field under the 

patch becomes : 

Cos (^) Cm (^) Gk(^) «(»>«(»> 

_ ^ q £ * (5) 

‘ ) 2^ 2* ( k 7 - k 7 ) 

QD m- 0 »t*0 \ ^ / 


The Neumann numbers e(m) and e(n) are , in general , represented as . 
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e(n) = [ 1 , |i = 0 ; 2 , p* 0 ] , and 1^= k 2 0 e r (l-j6) . The amplitude coefficient from (4) is : 


J _, ■ > 

< + ~ • ♦ J < *’ - *- > 


( 6 ) 


Here , the inner product is defined as : 

< > ■ // 

S 


(7) 


and < Y m , > = (ab) / ( e(m) e(n) ) 


2.2.2. Input impedance 

Since E z is assumed to be invariant to z , the impedance at the probe position is . 


d 





<*) 


where, I 0 is the probe current. However, the effect of the probe reactance, Xf , must be accounted for in 
order to obtain the overall input impedance at the probe position : 

X f c *o '1 (9) 

fr 

where t is the probe thickness and ti 0 is the free space impedance. 

Evaluating E z at x = x, , y = y p , from (5) gives 



( 10 ) 


.mixx 


ab 


coc ( )co« (- 


E E 

m-0 fbO 


fftTiX ft TtV W?ty 

' — )C06 (— —)CO* (- 

b 


— )€(m) €(«) 


i V - *- ] 


Therefore, 


DS 


~ J ^ d — £ £ 

dbt m m-0 ii-O 


mux. , 
cos ( )cos ( 


mux. 


— )cos (— — ^-)cos ( ^)e(w) €(«) 




( 11 ) 


[ “ - (! - J 'Qtot > w 1 


where, k 2 efr = e r ( 1 - 6 e ff)k 2 0 and 6 6 eff 


1/ Q TOT . Therefore, the input impedance seen at the probe is 


'*DS 


- Z 


DS 


jx f 


( 12 ) 


where, "DS" refers to double sum. 


2.3. Single-sum approach 
2.3.1. Electric field derivation 

Referring to figure 2, the patch is divided into two geometric regions: region I, where y>y p and region 

n, where y < y p . Solutions, for the electric field must be obtained for each region , in accordance with the 
inhomogeneous Helmholtz equation (equation (1) ). It is assumed that : 

£ * fx) g(y) ■ 2 f m (x) gjy ) ( 13 ) 

UNO 


where , f and g depend on the modal index m ( m - 0 , 1 , 2 ...) , and 




Figure 2 Two-Region , Single-Sum Patch Area 



with A,„ and B m being complex constants. The required boundary conditions are : 
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a x h 


m magnetic walls 


= o 


Enforcing the boundary conditions at x - 0 and at x - a gives : 

E, - £ A Cos (^) gjy) 

ukO a 


( 15 ) 


(16) 


Substitution of (16) into (1) , after multiplying each side of the new equation by Cos(mnx/a) , and 
integrating from x=0 to x=a , gives 

[ JL - . t > ] gjy) . ^ Ow(— £ ) » (y- y,) (17) 

dy* a a a 

The Green's functions solutions gjy) in the y-direction are continuous at the y = y p boundary. However, 
the derivatives at this boundary are not equal due to a jump condition from the probe excitation. From Felson 
et al. 1993, the components of the Green’s functions in regions I and II can be expressed as : 

gj y ) * A m fS y)fjy p ) , y>y, n & on 1 ( 18 ) 

gjy) - A m fj y ) fj y p ) , y<y f region n ( 19 ) 

where, f>(y) and f<(y) are functions belonging to regions I and II , respectively. 

Assuming sinusoidal variation , 


f y (y) = C m Cos (y y) * D m Sin (yy) , y*y p 


( 20 ) 
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/.(y) - E m Cos (yy) ♦ F m Sin (yy) , yzy p 


( 21 ) 


for C m , D m , E m and F m as complex constants. Using the Maxwell equation, which relates the curl of the 
electric field vector to the magnetic fields, gives : 


-ywn H z 


E ACos ( — ) 

m-0 & 


dgjy) 

dy 


( 22 ) 


The boundary conditions 




( 23 ) 


yield F m = 0 , and 


/< 00 - E m Cos (yy> 


(24) 


Letting A m ' = A ra E m , and realizing from (23) that C m = Cos(yb) and D m = Sin(yb) , we have : 


/> 00 * Coj[y(y-6)]) 


(25) 


Letting the Green's functions g,„(y) = G^tyj,,) , we have , for each respective patch region 

G ( y , y ) * A 'Cos( yy )Ow( y(y-i) ) , y*y , region 1 ^ 

y m P m p 

G ( y , y, ) - AjCos( yy )Gw( y (y p ~b )) , y±y f region U (27) 

Pm " 


Note that (26) and (27) are equal at the y = y p boundary of regions I and II . 
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Following the procedure outlined in Felson et al. 1993 to determine A,,, , we have 




Cos ( p ~) 

e(m) a 

a y Sin (yb) 


(28) 


Substituting (28) in (13), for the appropriate region , the single-index electric-field quantities of the microstnp 
antenna are : 


• VW Cos Cos Cm (Y>p ) Cos i yiy-y f )) . (2») 

'' ii aySm(.yb) a a 


- /<* H. 2, fW Co5 ( mnx ) cos (^£) Cos (YY) Cos (»(y,-»)) , »!f, W 

ii ay Sin (yb) a a 


The only unknown parameter that has not been dealt with so fitr is y - Y„ , the single-index propagation 
conam It is obtained from using either (26) or (27) for G^(y,y,) and substituting it into the homogeneous 

form of (1) : 


[ JL - (iH) 2 + k 2 ] A^Cos (yy)Cos (y(y-b)) * 0 (31) 

dy 2 a 


For non-vanishing cosines , 


N 


k 2 - (— ) J 


(32) 


where, k 2 = e r ( 1 - 6 eff )k 2 0 , as defined in section 2.2.1. 
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2.3.2 Input impedance 

Evaluating E z from (29) or (30) at the probe coordinate (x,, , y p ) and substituting into (8) , gives 




m nx 

_ c(w) co« : ( -) cos (l tj p ) COS (y m (y p - b)) 

a 


'ss 


m-0 


Y, sin (Y m b) 


(33) 


where, "SS" refers to single-sum Therefore, 

= ^ss * 


(34) 


where , X f is defined in equation (9). 

2.4. Comparison of the single-sum and double-sum input impedance expressions 

Mathematica 2.2 software was used to calculate all computer data. Sample antenna parameters chosen 
were a = 2" , b = 1.27 " , t = 20 mils , e r = 2.33, d = 62 mils , with a resonant frequency of 1 .93 GHz. The 
single sum input impedance for a probe position of x = 0.892”, y = b / 2" was simulated to six terms only, 
and held fixed in that manner. One of the double sum terms was also held fixed at six terms, and the 
remaining sum was separately computed by limiting the upper index to 5, 10, and 20 terms, respectively. This 
was performed to demonstrate the congruency of the single and double sum solutions. Figures 3a-c 
demonstrate this congruency on computer-generated Smith Charts. It is apparent that as the non-fixed 
double sum term approaches 20, the two solutions begin to overlap. Therefore, the single sum solution can 
serve as a replacement to the double sum solution, since it not only represents the double sum solution, but 
it also is easier to calculate. 

Table I is a computation of the execution times associated with the computer simulation of equations 
(12) and (34), for the antenna mentioned in the first paragraph of this section. The operating frequency 



2.33 , d • 62 mils ) 
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Convergence of Single and Double-Sum Input Impedance Loci 
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1 . 



(C) 


Figure 3 


Convergence of Single and Double-Sum Input Impedance Loci 


Table I - Input Impedance Execution Times for Double and Single Sum Terms 
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chosen was 1 .93 GHz. A reference to Table I shows that the double sum solutions take from 2 to 153 times 
as long to execute as compared to the single sum solutions. Again, the single sum solutions are desired for 
their lower computational times A graphical display of the congruency of double and single sum solutions 
is shown in figure 4a. This is further evidence that the two methods provide congruent results. The same 
simulation procedure mentioned above was performed for a second antenna with e r = 10.8, a = 0.9", b = 
0.6", d = 025", and a probe position of x = 0.54", y = b / 2", at f = 2 GHz. Results are shown in figure 4b, 

and are commensurate with those discussed above. 


3. Experimental verification of the data 

A Wiltron 360B Network Analyzer was used to test the rectangular-patch microstrip antenna for input 
impedance data, with respect to the S u input port parameter. The Smith Chart frequencies were swept from 
1.7 to 2.3 GHz. The test antenna was of patch dimensions a = 0.9", b= 0.6" , with e r =10.8, and substrate 
thickness of d = .025". The position of the probe was located at x = 0.54", y = b/2 ", relative to the bottom 
left comer (coordinate origin) of the patch. The same antenna dimensions and dielectric constant data were 
fed into the computer model to generate a Smith Chart. In the theoretical model, an operating frequency of 
2 GHz was used, which was the marked frequency on the Network Analyzer. 

Figures 5a and 5b are the theoretical and experimental versions, respectively. Usually, the main point of 
interest on either plot is the point which intersects the real axis, i.e., the point closest to the perfect 50C1 
probe-to-coax line matching point At 2 GHz, the unnormalized impedances of the input impedances at the 
probe location were (47.961 +j 1.891) Q and (47 842 +j 0.762) O for the theoretical and experimental 
data, respectively, when equation (9) was used to account for an approximation of the probe reactance. 
However, when equation (9) was omitted from the theoretical model, the theoretical probe impedance was 
equal to (47.961 + j 2.861), which is much better when compared to the experimental value. However, it 
should be noted that neglecting the probe analysis in theoretical modelling provides reasonably accurate 
values around the resonant frequency of the antenna, but it also moves the values of the theoretical input 
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Abs [Zin] 

(Double Sum and Single Sum Input Impedances) 



Frequency, GHz 


Figure 4 - Real Part of Double and Single Sum Input Impedances 
(a) - €,-2.33, d - .062", a- 2", b-1.27" 
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Theoretical Versus Eiperimental Data for Input Impedance Loci 
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impedances in a capacitive direction toward the lower half of the Smith Chart when frequencies are further 
away from the resonant frequency. The error associated with prediction of the imaginary part of the input 
impedance can be directly attributed to the fact that an approximation was used to forecast the value of the 
probe reactance. Only testing of the actual antenna provides the exact quantity. However, the use of 
equation (9) allows one to predict a more accurate value of the antenna input impedance with respect to 
predicting impedances over an entire given frequency span. The omission of equation (9) from the analysis 
has been found to cause the Smith Chart impedance data to be shifted downward. This means that the probe 
has a capacitive effect, which is not the actual case - its effect is an inductive one. 

Both Smith Chart curves are quite similar, and the small difference in results can be attributed to the fact that 
the theoretical model (cavity model) neglects consideration of surface waves and backlobe radiation (under 
the ground plane) due to reflection caused by an input impedance mismatch Furthermore, the theoretical 
analysis views the ground plane as being infinite in transverse (to z) directions. However, good results were 
still obtained, since it has been found (James and Hall 1981) that surface waves are negligible for thin- 
substrate patch antennas. Since the matching is reasonably close to 50 Q , reflection is minimal. Therefore, 
the analysis used was adequate enough to predict the experimental results. 

4. Conclusions 

The single-summation approach used to obtain the input impedance was shown to drastically reduce 
computational effort, when compared to similar double-summation calculations. The real value of this 
procedural economy is that the reduction in computation and cost does not sacrifice the accuracy needed 
Also, it should be emphasized that the double-summation approach required at least twice as much com - 
putation as that of the single-summation approach. When the double sum solutions are summed over 50 
terms each, it was found that it took 153 times as long to execute as that of an equivalent single summed 


term. 
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Abstract 

The need for high efficiency components has grown significantly due to the expanding role of fiber optic 
communications for various applications. Integrated optics is in a state of metamorphosis and there are many 
problems awaiting solutions. One of the main problems being the lack of a simple and efficient method of 
coupling single-mode fibers to thin-film devices for integrated optics. In this paper, optical coupling between a 
single-mode fiber and a uniform and tapered thin-film waveguide is theoretically modeled and analyzed. A novel 
tapered structure presented in this paper is shown to produce perfect match for power transfer. 

1. INTRODUCTION 

Integrated Optics has come a long way since early 1970’s. It is poised at the threshold of making a big 
impact in everyday world. Integrated optics has drawn into several disciplines, such as computers and 
microwave integrated circuit technology resulting in new fabrication technologies. The field is now in a state of 
flux and there are still many problems awaiting solutions. One of the main problems being the lack of a simple 
and efficient method of coupling optical fiber to thin-film devices for integrated optics. Although a number of 
connectors are available in the market, the efficiency and reproducibility are low. In this paper, a novel uniform 
and symmetrically tapered structure are analyzed mathematically and the results presented. 


There are certain papers published by researchers in the area of thin-film couplers that are worthy of note, in that 
they provide the background for the development of the proposed work and point towards issues raised by 
shortcomings of the previous work. Louisell [1] investigated broadband bi-directional couplers in which he showed 
that the phase constants and coupling coefficients vary with distance along two coupled transmission lines. 
Ulrich [2] has shown by analysis how light can be coupled into a thin-film by means of a prism-film coupler. Akira 
Ihaya [3] presented a mathematical model of a thin-film optical directional coupler consisting of a three-layered 
deposited glass films on the substrate, with coupling occurring between the first and the third film. Wilson and 
Teh[4] have shown a mathematical modeling of a tapered velocity directional coupler. Nelson [5] has theoretically 



examined the coupling of single-mode optical waveguides through the use of expanding and contracting tapers. 
Juichi Noda et al. [6] have shown a connection between single mode fiber coupler to a Ti diffused LiNo 3 strip 
waveguide. The authors Y.Cai et al.[7] have analyzed the coupling characteristics of a uniform structure. 


2. THEORETICAL BACKGROUND 

In this paper, we present a novel method of coupling light from a single mode fiber to two structures : Uniform 
and symmetrically tapered thin-film couplers. The model presented is an improved version of the uniform five 
layered structure presented by Y. Cai et al [7]. If a fiber is directly connected to the thin-film, as shown by Y.Cai et 
al [7]., it leads to a large mismatch of the field profiles at the interface and as a result, leads to loss of optical 
power. To overcome this optical reflection and radiation losses, they proposed a five layered structure which is 
shown in Figure L This consists of a coupling waveguide, buffer layers, thin-film layer and the substrate. To avoid 
mismatch of the fiber core and the coupling waveguide field profiles at the interface, the coupling waveguide 
dimensions are designed as 2aX2a, where 4 a’ is the radius of the fiber core. Guttmann et al[8]., have shown that 
the field distributions of the fundamental mode of the cylindrical fiber core and the garnet material waveguide 
differ slightly if the refractive index differences are smaller. Solgel solution is used at the interface such that any 
slight refractive index mismatch between the fiber core and the coupling waveguide is minimized. The buffer layer 
serves the purpose of coupling optical power to the thin-film waveguide of lower refractive index. Gadolinium 
Gallium Garnet material was used as the thin-film material because of the properties exhibited by the material. 
The buffer layers and the coupling waveguide are silica doped materials. The required refractive index of buffer 
and coupling waveguides is obtained by doping Ti0 2 with Si0 2 . 

The eigen mode equations of the five layered structure developed ( Figure 1 ) can be obtained by solving with 
w=g=2a, where ; w= the height of the coupling waveguide, a= radius of the fiber core, ngg=reffactive index of the 
coupling waveguide ( Garnet material ), n a =refractive index of the buffer layer, n w = refractive index of the thin- 
film material, k=wavelength of the wave ( 1.3pm ) and Wl= height of the thin-film coupling waveguide. The 
eigenmode equations developed by Y.Cai et al[9].,are : 


a = nlK^rxiK^.W / 2) ----- (1) 

Q, = y/(2x/Jl) 2 (£-n;)-Xi --- (2 ) 

K„{-l + y/tan 2 (g.X v ) + l} 

Q, = ( v \ ----- (3) 

tan(gA^, ) 

e, -n])-Kl ----- (4) 



The parameters k^. Q a and kgy, Q t are the x and y components of the wavenumbers. and Kg,, represent the 
sinusoidal variation in the x and y directions. The propagation constant of the coupling waveguide can be 
determined by solving equations ( 1 to 4 ), 


Pi = 6)2 Mo 


gx 


■K- 


(5) 


The propagation constant (3 W of the thin-film waveguide is obtained by replacing ilgg with n w , w with W1 and 
Pgg with p w in equations( lto5 ). The value of W1 for the thin-film is so chosen such that highest power 
coupling occurs between the coupling waveguide and the thin-film at this value. The optimum value W1 is chosen 
such that P ffi = Pw for no mismatch. The thickness of the buffer layer has a direct bearing on the power coupled 
from the coupling waveguide to the thin-film waveguide and an optimum value is chosen such that there is a 
maximum power transfer from the coupling waveguide to the thin-film waveguide. The propagation constant of 
the thin-film varies because of the change in value of w along the z direction. Therefore the mismatch, M depends 


on the propagation constant difference given by 


M(z) = 


P gg -Pw( z ) 

2C 


( 6 ) 


and the maximum fraction of power transferred between the coupling waveguide and the thin-film waveguide with 
a fixed mismatch M, is given by F 2 = ( 1 + M‘ ) '' ( 7 ) 

Application of the theory developed by Snyder et al[10] and Snyder [11], to the model analyzed in Figure 1, gives 
the power coupled into the thin-film for a uniform coupler as 

Z 

P’=PF 2 Sin 2 \{CIF)dz (8) 

0 

where P is the total power introduced, C the coupling coefficient between the center of the coupling waveguide 
and the center of the thin-film waveguide. In this paper, we present a novel uniform and symmetrically tapered 
structure ( Figure 2 ) where the modes of a uniform waveguide form a complete set and can propagate 
independently from one another, while the tapered modes are coupled together and adjust their characteristics to 
suit the varying transverse properties of the guiding structure as they are moving along the taper. During the 
process of coupling the light from the uniform coupling waveguide to the thin-film waveguide, the wave is kept in 
the lowest order mode. For the shape shown in this paper, ( Winn and Harris [12]), the power coupled at the 


narrow end is given by P 0 - P’ ( 1- F* Sin * J (C / F)dz ) ( 9 ) 

0 

In this paper while analyzing the theoretical model of the uniform and tapered coupler we made the assumption 
that the materials used are lossless and therefore there are no Fresnel’s reflections at the interfaces and that there 
are no lossess at the interface of the fiber core and the coupling waveguide. 



3. RESULTS AND DISCUSSION 


The material used for the analysis of uniform and tapered couplers is Gadolinium Gallium Garnet (GGG). The 
coupling waveguide material used is GGG. Buffer layer material is spun silica which is doped with Ti02. Thin- 
film material is polymerized solgel solution of Si0 2 and Ti0 2 . The refractive index of coupling waveguide, = 
1.9389 @ X = 1.3nm. The refractive index of the buffer layers is chosen, as n a = 1.9340 and the refractive index 
of the thin-film material, n w = 1.9450. Figure 3 shows the plot of C versus d, where C is the coupling coefficient 
and d is the distance from the center of the coupling waveguide to the center of the thin-film. This equation given 
by Snyder [11], is 


C 


FA 

2 7C. 


1/2 


ex pi 


-Vd 2 

4 P 2 


( 10 ) 


The refractive index of GGG versus wavelength is shown in Figure 4. The expression used to calculate the 


refractive index is 




( 11 ) 


where A t and Lj are the sellmeier coefficients,given by Wood and Nassau [13]. 

By solving the equations (lto5) with w=2a=g, the propagation constant (3gg = 9.3253804 pm' 1 was obtained. 

Figure 5 shows the propagation constant for the thin-film as a function of the film thickness Wl. The propagation 
constant of the thin-film equals the propagation constant of the coupling waveguide when the thickness 
W1=0. 886p.m. The maximum coupling power is attained at Pgg = Pw = 9.3253804 pm* 1 . Figure 6 shows the 
plot of coupled power versus the propagation distance for the uniform coupler of different d’s for perfectly 
matching constants p^ = Pw at thickness of thin-film w=0.886 pm. Figure 7 shows the power coupled for 
uniform and tapered couplers with respect to the distance z and the taper starts at length Ll= 700pm, with 
slopes (K) of the symmetrical taper at 0.0001, 0.00009 and 0.00002. Figure 8 shows the effect of the thickness of 
the thin-film and the role played by the taper in the power output. 


4. CONCLUSION 


The novel structure (Figure 2) exhibits broad-bandwidth coupling characteristics and it is easy to manufacture 
since it does not put serious constraints on the accuracy of the coupling length of the coupler. The taper has the 
distinctive advantage of confining the power within the taper such that the output stays approximately near the 
value of the power introduced at the start of the taper. The outstanding feature of the taper is the higher efficiency 
as compared to the more conventional devices. 




Figure 1. Schematic diagram of Single-mode fiber to 
thin -film waveguide using a coupling waveguide 
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Figure 3. Coupling coefficient versus sepe ration distance 
between center of coupling waveguide to center of thin- 
film waveguide. 



Figure 2. Schematic diagram of the proposed novel 
tapered structure. 
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Figure 4. Refractive index of gallium garnet as a function 
of wavelength. 



FigureS. Propagation constants versus thickness of thin- 

film. 


Figure 6. Plot of power coupled into thin-film versus 
length for various values of d. 
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Figure 7. Comparison of uniform power and tapered 
power coupled Into thin-Glm with taper at LI *700 jun fo r 
different cases of K for a thin-film thickness of 
Wl- 0.886 jam. 



Figure 8. Plot of coupled power versus propagation 
distance for a taper stucture wtth Ll-1700 ^im and a 
slope K~0.00025 for different thickness of thin- film- 
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